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AU3TRAC1 


The  vertical  accelerometers  were  used  to  make  a  six  month  (January  to  July  1971)  data 
collection  run  which  is  being  used  to  determine  detection  thresholds  and  azimuthal 
dependences.  This  run  was  made  at  a  surface  site  at  Camp  Elliott,  San  Diego.  The 
horizontal  accelerometers  were  completed  and  placed  in  operation  and  data  gathered 
simultaneously  with  the  vertical  accelerometers  for  several  months  at  a  surface  site 
at  Camp  Elliott,  San  Diego.  Both  a ( vertical  and  horizontal  were  used  in  a  mine  test 
at  5300  and  6700  feet.  The  results  showed  that  the  major  improvement  came  in  the 
horizontal  data  which  showed  no  diurnal  effect  as  is  seen  in  the  surface  data.  The 
absolute  values  of  the  horizontal  mine  data  taken  at  5300  feet  both  for  day  and  night 
in  the  seismic  band  (45  second  period  and  shorter)  were  as  good  as  the  quiet  night 
surface  values  at  San  Diego.  There  was  not  much  difference  in  the  vertical  data  taken 
at  depth  and  at  the  surface. 
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SUMMARY 


The  principal  results  of  this  contract  are  the  followings 

1.  The  vertical  accelerometers  were  used  to  make  a  six  month  (January  to 
July,  1971)  data  collection  run  which  is  being  uaed  to  determine  detec¬ 
tion  thresholds  and  azimuthal  dencndences.  This  run  was  made  at  a  sur¬ 
face  site  at  Camp  Elliott,  San  Diego. 

2.  The  horizontal  accelerometers  were  completed  and  placed  in  operation  and 
data  gathered  simultaneously  with  the  vertical  accelerometers  for  several 

months  at  a  surface  site  at  Camp  Elliott,  San  Diego. 

* 

3.  Both  a  vertical  and  horizontal  were  used  in  a  mine  test  at  5300  and  6700 
feet.  The  results  showed  that  the  major  improvement  came  in  the  horizon¬ 
tal  data  which  showed  no  diurnal  effect  as  is  seen  in  the  surface  data. 
The  absolute  values  of  the  horizontal  mine  data  taken  at  5300  feet  both 
for  day  and  night  in  the  seismic  band  (45  second  period  and  shorter) 
were  as  good  as  the  quiet  night  surface  values  at  San  Diego.  There  was 
not  much  difference  in  the  vertical  data  taken  at  depth  and  at  the  sur¬ 
face  . 
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Introduction 


The  Final  report  on  ARPA  Contract  IGPP/AF/F-4M620-71-C-0128  is  divided 
into  four  sections.  The  first  two  will  deal  primarily  with  data  taken  at 
the  Star  Mine  at  Wallace,  Idaho  and  at  Camp  Elliott,  San  Diego,  California. 
The  last  two  sections  deal  primarily  with  instrumental  improvements  in  the 
vertical  accelerometer  and  a  first  description  of  the  horizontal  accelero- 
meter.  These  last  two  sections  have  been  submitted  for  publication. 

Section  I.  Operation  of  Vertical  and  Horizontal  Accelerometers 
at  Croat  Depth. 

Section  II.  Teleseismic  Signals  Eetween  and  Below  the  Microseism 
Peaks.,  Bartholomew,  B.,  Block,  B.  and  Dratler,  J.,  Jr, 

Section  III.  Improvements  in  the  Wide  Band  Vertical  Quartz 
Accelerometer.  Block,  B.  and  Dratler,  J.,  Jr. 

Section  IV.  A  Wide  Band  Horizontal  Accelerometer  with  Preliminary 
Earth  Normal  Mode  and  Seismic  Investigations.  Dratler,  J.,  Jr., 
and  Block,  B.  (In  press,  Geophys .  Jour .  Roy .  Astro .  Soc  . ) 
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Section  I.  Operation  of  vertical  and  horizontal  accelerometers  at 
great  depth. 

Data  was  obtained  in  the  Star  Mine,  Wallace,  Idaho  at  depths  of  5300  and 
6700  feet  for  both  vertical  and  horizontal  broad  band  quartz  torsional  accele¬ 
rometers.  The  experimental  set-up  will  bo  described  for  the  station  at  5300 
feet  as  it  was  typical. 

The  Star  Mine  is  a  zinc,  lead  and  silver  producing  mine.  The  matrix  rock 
is  quartzite.  The  location  of  the  seismic  observatory  was  chosen  to  be  as 
far  as  possible  from  the  active  excavation  sites  (unfortunately  only  a  few 
hundred  feet)  and  to  be  located  in  a  tunnel  which  was  not  used  in  the  mine 
ventilation  system.  The  tunnel  was  closed  off  by  means  of  metal  doors  at 
either  end.  At  5300  feet  the  temperature  in  the  tunnel  was  95°  F  with  100% 
humidity  present.  The  tunnel  showed  relatively  little  rock  damage  compared 
to  other  possible  sites  investigated  but  every  day  in  the  1000  feet  from  the 
shaft  elevator  to  the  instruments  several  small  rockfalls  were  found.  This 
mine  has  a  history  of  rock  bursts.  In  some  locations  the  rock  damage  was  so 
extensive  that  the  walls  had  caved  in  two  feet  in  a  height  of  10  feet  even 
when  held  back  by  steel  sheets  and  rock  pins.  The  acceleration  sensors  were 
housed  in  a  timber  shack  built  in  a  bedbug.  The  shack  was  covered  with  bur¬ 
lap  backed  by  polyethylene  film.  No  temperature  control  was  applied  in  the 
shack.  The  shack  was  used  merely  to  keep  drafts  off  the  instrument  and  to 
protect  against  rockfalls. 
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A  pier  was  installed  by  chipping  away  the  shattered  base  rock  and  using 
a  thin  layer  of  grout  to  lay  a  granite  slab.  A  100  feet  of  cable  was  installed 
to  another  bedbug  where  a  similar  shack  was  built  to  house  the  electronics  and 
data  logging  equipment.  This  shack  was  cooled  by  an  air  conditioner  and  a  de¬ 
humidifier  operated  continuously.  This  cooling  and  dehumidiiication  was  found 
necessary  to  keep  the  electronics  and  the  data  logger  operating.  A  few  words 
about  the  data  logging  system  are  in  order.  This  is  a  twenty  channel,  9  track 
800  BPI  system  made  by  Monitor  Laboratories,  San  Diego.  The  system  operated 
successfully  except  for  a  power  supply  failure  probably  induced  by  the  shack 
cooling  system  failing.  Continued  operation  occurred  after  the  mechanically 
traumatic  move  into  the  mine  site  at  6700  and  5300  feet.  The  tape  deck  stayed 
within  specifications  at  all  times  and  the  digital  electronics  operated  success¬ 
fully.  Loss  of  components  occurred  in  the  accelerometer  electronics  when  the 
shack  cooling  and  dehumidifying  system  failed.  As  soon  as  these  systems  were 
re-established  giving  a  temperature  of  80°  F  and  a  humidity  level  somewhat 
.less  than  100°  the  electronics  problems  ceased. 

Both  digital  and  analog  records  were  taken  over  several  weekends  in  the 
tidal,  earth  normal  mode  and  seismic  filter  channels  These  channels  are  de¬ 
scribed  in  detail  in  Section  III  and  IV,  The  sampling  rate  was  once  a  second 

4 

with  a  digitizer  dynamic  range  of  10  in  each  filter  channel. 

Data  in  this  mine  can  be  taken  only  over  the  weekend  from  Saturday  morning 
to  Sunday  night  because  the  mine  is  on  continuous  three-shift  operation  through¬ 
out  the  remainder  of  the  week.  The  data  for  the  5300  foot  station  will  be  dis¬ 
cussed  below.  The  data  from  the  6700  foot  station  was  much  noisier  than  that 
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obtained  at  5300  feet.  This  is  attributed  to  the  massive  rock  damage  and  con- 

# 

sequent  motion  at  this  depth  in  the  mine.  The  walls  were  literally  caving  in 
at  6700  feet.  This  field  operation  could  not  have  been  carried  out  without  the 
complete  cooperation  of  personnel  at  the  Star  Mine  site  and  the  management  of 
the  Hecla  Mining  Company. 

0 

Figure  1  shows  typical  horizontal  data  obtained  from  the  5300  foot  level 
of  the  Star  Mine  taken  on  a  weekend.  The  seismic  filter  response  functions  of 
both  a  horizontal  and  vertical  accelerometer  are  found  in  Sections  (II-IV). 

For  comparison,  typical  data  taken  at  a  surface  site  at  Camp  Elliott,  San  Diego 
are  also  given  in  Fig.  1.  In  essence,  what  Fig.  1  shows  is  that  the  vertical 
absolute  noise  at  both  sites  is  about  the  same  through  this  seismic  filter. 

The  horizontal  noise  at  Camp  Elliott  varies  significantly  from  day  to  nigh:. 

The  daytime  noise  level  is  almost  always  being  much  larger  than  the  night¬ 
time  noise  level.  However,  the  horizontal  noise  at  the  5300  foot  level  does 
not  show  this  diurnal  change  but  remains  appreciably  constant  throughout  a 
weekend.  The  absolute  noise  level  at  5300  is  about  the  same  s.ize  as  an 
average  day  at  Camp  Elliott.  But  what  happens  at  seismic  frequencies  (i.e., 
defined  as  a  high  pass  filter  with  6  db  point  at  45  seconds)  is  seen  to  be 
that  the  seismic  horizontal  noise  level  at  5300  feet  is  always  about  the 
amplitude  of  the  best  nighttime  seismic  horizontal  noise  level  at  Camp  Elliott. 
For  the  daytime  data,  the  5300  foot  horizontal  noise  level  is  considerably 
smaller  than  the  daytime  Camp  Elliott  data.  (See  Fig.  1  Caption)  Two  typical 
earthquakes  taken  at  5300  feet  are  shown  in  Fig,  2  for  the  vertical  accelero¬ 
meter  along  with  their  high  pass  filtered  data.  Information  on  these  earth¬ 
quakes  is  contained  in  the  caption  of  Fig.  2.  This  data  was  hand  digitized 

-  5  - 
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from  analog  records  when  the  data  logger  was  down.  It  can  be  done  to  some 
extent  but  one  experience  convinces  the  author  that  this  way  of  proceding 
is  hard,  costly  and  loses  much  dynamic  range. 
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Figure  1  Caption 


Data  from  the  horizontal  Instrument  showing  day  end  night  variations 
at  the  Camp  Elliott  Station  (ground  level)  and  the  absence  of  such 
variations  at  Wallace,  Idaho  (5300  feet  below  the  surface).  The  data 
are  also  shown  after  being  convoluted  with  a  hipass  filter  that  cuts 
af  45  sec. 

a)  Raw  data  from  Wallace:  S.T.  13:29:10  local  time. 

b)  Ras  data  from  Elliott:  S  T.  12:04:08  Local  time. 

c)  Filtered  data  from  Wallace:  S.T.  13:29:10  local  time. 

d)  Filtered  data  from  Elliot:  S.T.  12:04:08  local  time. 

e)  Filtered  data  from  Elliott:  S.T.  00:04:08  local  time. 

f)  Filtered  data  from  Wallace:  S.T.  01:29:10  local  time. 

g)  Raw  data  from  Elliott:  S.T.  00:04:08  local  time. 

h)  Raw  data  from  Wallace:  S.T.  01:29:10  local  time. 
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figure  2  Caption 


Two  earthquakes  seen  through  the  vertical  seismic  filter 
recorded  on  a  stripchart  at  the  5300  foot  level  of  the  star 
mine  are  shown  in  A  and  C.  The  result  of  hand  digitizing 

*nd  subsequent  high  pass  digital  filtering  (6  db 
point  at  45  seconds)  is  seen  in  B.  and  0. 

• 

Event 

Location 

Origin  Time  (GMT) 

Magnitude  ^B) 

Distance 

A 

North  Atlan¬ 
tic  Ridge 

17.05 

9/4/71 

4.6 

66* 

C 

5.3#S,  153.0E 
New  Ireland 

04:35:35 

9/5/71 

4.5 

95* 
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- - A-** - |  pVt'1  • ' 


-aM  fMl  .  .,'!.!!'<•  •»•  ,1 


B 


JL 

2  VOLTS 


,1  v. 


T 


VV'^-‘V — AV- 


— >  10  MIN 


h- 


-  10  - 


11 


Section  II. 


TELESEISMIC  DETECTION  BETWEEN  AND  BELOW  THE  MICROSEISM  PEAKS 


Teleseisnic  signals  and  background  noise  will  be  examined  in 
the  frequency  band  .02  -  .1  Hz  (10-50  seconds)  extending  from  the  noise 
valley  near  25  seconds  to  the  noise  minimum  or  crevice  between  the  micro¬ 
seism  peaks.  Power  spectra  and  digital  band  pass  filtering  of  the  back¬ 
ground  noise  and  two  small  events  will  be  given  by  way  of  illustration. 

The  data  from  the  broad  band  sensor  (the  Block-Moore  quartz  fibre 
12  3 

accelerometer,  *  *  Camp  Elliott,  San  Diego)  are  passed  through  an  active 
seismic  prefilter  (the  total  transfer  function  of  the  system  is  shown  in 
Fig.  1)  and  digitally  recorded  on  magnetic  tape,  yielding  a  dynamic  range 

4 

of  10  with  the  seismic  prefilter  bandwidth.  Other  means  of  recording  the 
data  can  cause  a  severe  degradation  of  either  dynamic  range,  band  width,  or 
both,  leading  to  impaired  detection  capability.  For  instance,  analogue 
recording  (e.g.  photographic)  of  seismic  data  using  a  notch  filter  to  reduce 
the  visual  contribution  of  the  large  6  second  microseisms  can  degrade  the 
dynamic  range  by  a  factor  of  100  from  this  recording  method  alone.  The 
effect  of  the  wings  of  the  notch  filter  in  lowering  the  signal  level  below 
the  recorder  threshold  will  further  degrade  the  data.  The  usefulness  of 
the  data  between  the  microseism  peaks,  where  both  explosions  and  earthquakes 
possess  high  signal  levels,  is  thus  severely  limited.  Subsequent  digitizing 
of  such  data  cannot  recover  this  loss  in  dynamic  range  and  bandwidth. 

The  event  in  Figs.  2  and  3  was  recorded  simultaneously  by  a  vertical 
and  .horizontal  instrument  in  the  surface  vault  at  Camp  Elliott.  The  d.vta 
shown  in  Figs.  2  and  3,  seen  through  the  response  function  of  Fig.  1,  are 
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from  an  event  near  Oaxaca,  Mexico  (N  15.1°,  W  96.3°)  of  magnitude 
Mb  =  4.7  on  June  6,  1971  02:10:26.9  GMT  at  a  depth  of  33  km.  The  angular 
distance  A  from  the  sensor  was  25.9°.  The  angle  a  between  the  great 
circle  path  connecting  the  epicenter  with  Camp  Elliott  and  the  axis  of 
maximum  sensitivity  of  the  horizontal  instrument  was  83°, 

o 

In  these  figures  the  raw  data,  denoted  by  A,  can  be  compared 
with  the  narrow  band  passed  signals  at  10  seconds  (.05  Hz  filter  width), 
15,  20,  and  30  seconds  (.01  Hz  filter  width),  denoted  by  B,  C,  D,  and  E 
respectively.  The  power  spectra  of  the  event  (dotted  line)  and  background 
noise  (solid  line),  computed  from  20  minute  intervals  during  and  just  be¬ 
fore  the  earthquake,  are  shown  in  Figs.  4,  5.  On  this  day,  the  18  second 
microseism  peak  for  vertical  acceleration  is  low  and  the  noise  level  in 
the  crevice  approaches  the  level  of  the  noise  valley  between  20  -  50  se¬ 
conds,  presuming  the  least  count  noise  of  the  digitizer  is  below  the. 
background  level,  a  condition  easily  met  with  modern  digitizers.  The 
horizontal  data  (on  Fig.  5)  show  a  similar  maximum  in  the.  signal  to 
noise  ratio  near  15  seconds,  but  the  signal  degrades  rapidly  until  there 
is  no  discernable  signal  at  30  seconds. 

The  noise  crevice  may  vary  in  amplitude  by  a  factor  of  3  above 
the  level  shown  in  Fig.  4,  although  for  a  majority  of  the  time  it  is 
clearly  discernable,  and  the  best  place  to  take  narrow  band  data. 

The  raw  data  from  the  vertical  accelerometer  for  a  second 
small  earthquake  near  Oaxaca,  Mexico  (N  15.2°,  W  94.6°)  of  magnitude 


Mb  =  5.2  on  February  20,  1971  12:25:01.7  GMT  appears  In  Fig.  6,  The 

earthquake  was  71  kt«.  deep,  the  angular  distance  to  the  sensor  i  a  27.0“ 

and  a  =  80°.  The  power  spectra  for  this  event  (Fig.  7  )  show  that  the 

crevice  noise  level  is  about  7  db  above  the  crevice  level  in  the  previous 

example,  but  the  signal  to  noise  ratio  between  12  -  15  seconds  is  still 

good. 

The  advantaCes  of  digital  collection  and  processing  of  data 
compared  Co  the  irreversible  losses  inherent  in  analogue  collection  of 
data,  even  with  subsequent  digitization,  can  make  the  difference  between 
detecting  and  missing  an  event. 


Figure  1.  The  power  spectra  of  the  response  functions  for  the  vertical 
and  horizontal  accelerometers  through  the  active  seismic  prefilter.  The 
response  functi  ns  in  (Volts/Micron)2  are  plotted  in  db  (0  db  =  KVolt/Micron) 
as  a  function  of  the  period  in  seconds.  The  microns  refer  to  microns  of 
ground  motion. 

Figure  2.  Data  from  the  vertical  instrument  for  the  event  on  June  6,  197i, 

The  ruw  digitized  output  from  the  accelerometer  and  active  seismic  pre- 
filter  is  denoted  by  A,  The  narrow  band  filtered  outputs  at  10,  15,  20, 
and  30  seconds  are  denoted  by  B,  C,  D,  and  E  respectively.  The  filter 
width  is  .05  Hz  at  10  seconds  and  .01  Hz  at  15,  20,  and  30  seconds.  The 
start  time  for  all  scries  is  02:1M:49  GMT  June  6,  1971, 

Figure  3.  Data  from  the  horizontal  instrument  for  the  event  on  June  6,  1971. 
The  raw  digitized  output  from  the  accelerometer  and  active  seismic  pre- 
filter  is  denoted  by  A.  The  narrow  band  filtered  outputs  at  10.  15,  20, 
and  30  seconds  are  denoted  by  B,  C,  D,  and  E  respectively.  The  filter 
width  is  .05  Hr.  at  10  seconds  and  .01  Hz  at  15,  20,  and  30  seconds.  The 
start  time  for  all  scries  is  02 : iu -t»9  CHT  June  6,  1971. 

Tigure  U.  Power  spectra  of  the  raw  vertical  data  from  Fig.  2  for  the 
earthquake  on  Juno  6,  1971.  The  vertical  scale  is  in  db  (o  db  *  1  aicron2/hzJ 
The  microns  refer  to  microns  of  ground  motion.  The  dotted  lino  was  cocputcd 
from  a  20  minute  interval  during  the  earthquake,  and  the  solid  line  is  from 
a  20  minute  interval  just  before  the  event. 
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Figure  5.  Power  spectra  of  the  raw  horizontal  data  from  Fig.  3  for  the 

•  2 

earthquake  on  June  6,  1971.  The  vertical  scale  is  in  db  (o  dl  -  1  micron  /Hz) 
The  microns  refer  vo  microns  of  ground  motion.  The  dotted  line  was  computed 
from  a  20  minute  interval  during  the  earthquake ,  and  the  solid  line  is  from 
a  20  minute  interval  just  before  the  e^ent. 

Figure  6.  Output  from  the  vertical  accelerometer  on  February  20,  1971 
showing  the  event  near  Oaxaca,  Mexico,  and  the  section  of  noise  used  to 
calculate  the  background  noise  power  spectrum. 

Figure  7.  Power  spectra  for  the  raw  vertical  data  from  Tigure  6  for  the 
earthquake  on  February  20,  1971.  The  vertical  scale  is  in  db  (o  db  * 

1  aicron^/hz).  T^c  dotted  line  was  computed  fr^n  a  3C  minute  interval  dur¬ 
ing  the  earthquake,  and  the  solid  lino  is  from  a  30  minute  interval  just 


before  the  event. 
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IMPROVEMENTS  IN  THE  WIDE  BAND  VERTICAL 
TORSION  ACCELEROMETER 
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Institute  of  Geophysics  and  Planetary  Physics 
University  of  California,  San  Diego 


IMPROVEMENTS  IN  THE  HIDE  BAND  VERTICAL  TORSION 
ACCELEROMETER 


Significant  improvements  have  been  made  in  the  wide  band 

vertical  torsion  accelerometer  since  the  original  design  and 

1  2 

operation  were  reported  upon.  * 

In  this  paper,  these  improvements  are  outlined  and  some 
typical  data  taken  with  the  modified  instrument  are  shown. 

Instrumental  Improvements 

The  vertical  accelerometer  is  compesed  of  a  horizontally 
stretched  fiber  in  torsion  which  supports  a  mass  on  a  moment  arm 
against  gravity.  The  mass  in  the  shape  of  a  flat  plate  forms 
the  center  element  in  a  balanced  parallel  plate  capacitance 
position  sensor.  The  position  sensor  uses  a  phase  sensitive 
detector.3  The  mechanical  Q  of  the  instrument  is  high  so  as 
to  lower  the  Brownian  noise  in  the  bandwidth  of  interest,  since 
the  mass  used  (approximately  10  grams)  is  small. 

The  instrument  is  small  (8  inches  OD)  and  comes  with  its 
own  pressure  and  temperature  controlled  environment.  It  is 
suitable  for  borehole  use. 

Several  modifications  were  made  in  the  capacitance  position 
sensor  to  improve  the  signal-to-noise  ratio  directly.  The  capa- 
citor  plate  area  (Fig.  1)  was  increased  from  9.5cm  to  19.5cm 
and  the  spacing  from  the  position  detector  mass  plate  to  either 
capacitor  plate  was  decreased  from  1.5mm  to  0.5mm.  These  changes 
result  in  an  overall  increase  in  signal  level  of  6.1  without  change 
in  any  noise  level.  The  free  period  has  been  lengthened  to  a 


typical  value  of  0.8  second  from  the  previous  value  of  about  0.6 
second.  This  was  made  possible  by  increasing  the  fiber  pulling 
speed  which  allows  a  more  uniform  diameter  fiber  to  be  pulled 
from  the  matrix  rod  while  keeping  the  overall  length  constant. 

In  the  new  models,  the  drive  voltage  has  been  reduced  from  40 
volts  peak  to  peak  at  the  capacitor  plates  to  15  volts  peak  to  peak. 
This  reduction  in  drive  voltage  diminishes  the  non-linear  destabilizing 
force  by  a  factor  of  7.  The  increase  in  plate  area  and  the  decrease 
in  plate  separation  together  increase  the  destabilizing  force  by  about 
a  factor  of  18.  The  net  increase  in  destabilizing  force  is  then  a 
factor  of  2.5,  which  has  not  been  found  troublesome.  The  net  increase 
in  signal  level  from  the  above  sources  is  a  factor  of  2.3. 

The  capacity-to-ground  CDQ  of  the  position  detector  plate 
(Fig.  2,  D1  and  Fig.  3)  was  reduced  by  modifying  the  center  plate 
clamping  assembly  (Fig.  2,  D4)  with  an  insulating  quartz  flat 
(Fig.  2,  D2)  upon  which  the  detector  plate  is  clamped.  A  1  mill 
tungsten  wire  is  used  as  the  signal  lead  directly  from  the  detec¬ 
tor  plate.  Another  1  mill  tungsten  wire  is  used  to  ground  the 
remainder  of  the  center  plate  clamping  assembly.  The  capacity 
CDG  *n  the  neM  vers*on  of  t*ie  instrument  would  have  been  large 
compared  to  that  in  the  published  version^ ^  because  of  the  intro¬ 
duction  of  eddy  current  damping  ears  (Fig.  2,  D3).  The  use  of 
these  damping  ears  in  controlling  the  Q's  of  unwanted  modes  is 
discussed  below.  As  a  general  practice,  the  reduction  of  C_„  from 
whatever  source  is  advisable  as  CDQ  acts  as  a  voltage  divider  to 
reduce  the  signal  level  at  the  preamplifier.  To  this  end,  the 
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preamplifier  (PAR  Hodel  11?)  i *  built  so  that  its  input  ITT  Is  bunted 

Sirectly  onto  the  seal-through  in  the  outer  vacuum  can  of  the  instrument. 

The  protective  dic^  j  at  the  input  ITT  are  removed  to  reduce  C  . 

DC* 

Modifications  suggested  by  the  manufacturer  were  made  to  the  mixer  of  the 
phase  sensitive  detector  (PAR  Model  120)  to  reduce  its  noise  level. 

>r  damping  certain  modes  of  the  maas-torsion  fiber  system*  it 
was  found  that  the  damping  ears  (Tig.  26  D3)  were  necessary  The 
snail  magnets  (rig.  4)  provide  a  field  which  is  perpendicular  to  the 
pl.n.  of  th.  damping  aars.  Thi.  configuration  tap.  out  th.  tvi.ting  rod.  -ho,, 
axis  is  along  the  ana  which  supports  the  mass  and  also  the  translational 
nodes  which  are  perpendicular  to  the  mass  arm  axis.  The  Q's 
or  the  other  modes  are  damped  by  the  magnets  mounted  with  oppo¬ 
site  poles  facing  on  the  two  detector  plates  Ir  prac- 

Uc,  the  Q  s  of  the  nodes  are  determined  by  the  magnetic  fields 
present  and  not  by  the  residual  atmosphere  in  the  inner  high 

acui'm  can.  The  Q  of  the  fundamental  mode  was  varied  from  20 
to  200  in  various  models.  No  effects  in  the  earth  nomal  mode  or 

seismic  channels  described  below  were  seen  coming  fren  instrument  code  excita¬ 
tion  and  nonlinear  coupling  to  the  fundamental  when  the  ambi^t 
acceleration  field  was  produced  by  nearby  jetc  it  Miramar  vaval  Air 

ftatioa  and  could  even  be  felt  causing  ground 'Vibration. 

All  stainless  steel  parts  in  the  center  clamp  assembly  ver* 

replaced  by  beryllium  copper  parts.  This  was  to  ensure  that  no 

effect  would  cone  from  local  changing  magnetic  fields  cr  magnetic 

storms.  The  mass  was  made  of  OFNC  copper.  A  magnetic  field  wi«h 
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a  60  second  period  having  an  amplitude  of  103  game  at  the  accel¬ 
erometer  produced  no  visible  effect  on  the  seismic  channel  which 

has  the  highest  gain  at  the  60  second  period  of  any  output  channel. 

The  improvements  in  the  temperature  control  system  involve 

both  the  passive  and  active  parts  of  the  system.  The  accelero¬ 
meter  (Fig.  5)  is  housed  in  a  hard  vacuum  (*10"7mHg)  stainless  steel 
can  sealed  with  a  gold  0-ring.  This  in  turn  In  housed  in  an 
aluminum  vacuum  can  (<l«Hg)  sealed  with  a  rubber  0-ring.  This 
outer  vacuum  can  is  hung  from  a  rigid  thermal  insulating  structure 
made  from  fiber  glass  (Fig.  5).  An  oute?  aluminum  surf ace. called 
the  control  surface .end  open  on  the  top  supports  the  fiber  glass 
structure.  This  control  surface  has  the  heater  winding  ar.d  sensing 
thermistor  (Figs.  5  and  6)  on  its  outer  surface.  The  control  eur 
face  along  with  the  fiber  glass  structure  forms  an  additional  thermal 
filter  from  the  outside  world.  The  control  surface  sits  upon 
legs  made  of  thin  wall  stainless  steel  tubing  for  thermal  iso¬ 
lation  from  the  tilting  platform.  The  control  surface  is  nested 
inside  a  container  and  all  spaces  are  filled  with  polystyrene 
balls.  In  practice,  the  tilting  platform  Is  covered  with  Innulc- 


An  alternate  heat  path  to  ,^hat  through  the  nested  series  of 
thermal  filters  described  above  is  through  the  electrical  cables 

to  the  accelerometer  and  temperature  sensor. 

In  order  to  maximise  the  effective  thermal  resistance  of 

these  cables,  lengths  of  them  are  coiled  inside  the  thermal 
enclosure  before  they  emerge  into  the  outside  world.  In  addition, 
lightweight  cables  are  used  to  connect  the  mechanical  system  to 
the  external  electronics.  The  coaxial  cables  are  teflon-insulated 
FG187/U,  while  the  other  conductors  are  part  of  a  30AWG  ribbon 
cable.  After  emerging  from  the  thermal  enclosure,  these  cables 
are  connected  to  heavier  standard  cables  (RG2/U  and  18AWG  multi¬ 
conductor),  which  connect  the  deployable  unit  to  the  electronics 
par.el. 

With  the  thermal  insulation  constructed  as  outlined  above, 
the  time  constant  for  heat  transfer  between  the  outside  world  and 
the  fiber  is  several  hours.  The  thermal  filtering  drastically 
reduces  the  effects  of  external  temperature  fluctuations  with 
periods  of  a  few  minutes  or  less.  However,  the  high  DC  gain  of 
the  temperature  controller  eliminates  the  low  frequency  temperature 
variations. 


5 


The  fed back  temperature  contoller  is  a  proportional  control 
system  with  extremely  high  DC  gain.  A  thennistor  measures  the  tespet ature 
at  a  point  on  the  control  surface  (Tig.  5).  The  temperature  sensor  bridge 
and  preamplifier  compare  the  thermistor  resistance  to  that  of  a  standard. 
An  error  signal  is  produced  proportional  to  their  difference. 

This  error  signal  is  sent  to  the  fedback  power  amplifier 
whose  output  determines  the  current  in  the  resistive  heater 
wound  around  the  control  surface.  The  error  signal  changes  the 

heater  current  to  bring  the  control  surface  temperature  beck  to 
the  null  point.  Since  the  null  point,  which  is  determined  by 

the  resistance  of  the  standard  in  the  ssnaor,  exceeds  room 

temperature,  there  is  always  current  in  the  heater  winding.  This 
\ 

mean  curiont  compensates  for  thermal  leakage  through  the  inr* la- 
tion.  Fluctuations  in  the  heater  current  maintain  the  texpera- 
ture  at  the  thermistor  fixed  within  »»  x  X<T*»C  per  100  sec 
and  «  x  10"8oC  per  day.  Theee  limits  are  obtained  by  assuming 
the  output  in  the  seismic  filter  rmd  the  drift  In  the  tidal  chan¬ 
nel  are  caused  ty  temperature  fluctuations  alone. 

All  models  of  these  instruments  have  extremely  low  drift. 

The  first  one  was  operated  continuously  for  a  year  and  has  shown  it-s 

than  10  *g  drift  during  this  time.  A  conservative  upper  bound  to 

-10 

the  drift  rate  is  thus  10  g/day. 

In  practice,  the  active  temperature  controller  and  the  passive 
thermal  insulation  interact  considerably.  Hot  only  does  the  in* 
sulation  determine  the  ncan  heater  current  needed  to  maintain  a 
given  temperature,  hut  it  also  affects  the  stability  of  the  temper »- 
turc  controller.  Tor  small  fluctuations  in  the  heater  currr.it  ,  u.«-  in- 
-.ulvion  hf’lps  to  control  the  hea*  loss  from  the  control  surface  »n?  •  ,u.. 
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the  thermal  lag  between  the  heater  and  themie'.or.  This 
the real  lag  is  primarily  responsible  for  difficulties  in  stabili 
sing  recalcitrant  control  system.  As  in  the  case  of  the  design 
3f  fhe  lnsul**;ion  itself,  calculations  of  the  stabilising  effects 
nre  rade  prohibitively  difficult  by  the  fact  that  the  parameters 

are  distributed.  Such  improvements  are  moat  quickly  made  by  trial 
and  error. 

tb,  ta^.raturo  controller  can  Imthh  th.  Mod  for  tkorMl 
insuLtlon  by  providing  .noil  thoro.1  fluctuation,  of  lto  » 

Tie  thermal  insulation  between  the  control  surface  and  the  outer 
aluminui:  can  forestalls  this  problem.  Being  a  lon-linear  system, 
the  thermal  controller  is  subject  to  marginal  instabilities. 

Usually t  th...  «  r.th.r  high  froguancioo  wh.ro  th.  controll.r 

loup  gain  1.  low.  th.  prosenc.  of  strong  onvlronnontnl  tempera 
tur.  drift,  howavar,  th.  controll.r  .om,1m.  cycl„  „  tll. 

higher  tol.Ml.nIc  froquaoci...  To  uduu  the  off.et  of  the.. 

low-l.v.1  thermal  osclllatlonc,  the  control  .urf.c.  1.  fl.c.d 

on  tll.  outside  Of  th.  soft  vacuum  can  and  Is  ito.lf  insulated 
fro-  direct  external  te.por.tur.  varlotlons.  The  heat  c.paclty 
of  tho  control  surface  Is  significantly  lower  then  that  of  ,h, 
soft  vacua  can.  thus  .taking  the  te.pen.ture  controll.r  easier 
to  stabilize.  The  pres.nt  isolated  control  surface  is  open  a, 
top.  At  low  frequencies,  themal  gradient,  can  exist  between 
the  control  surface  and  the  soft  vacua  can.  Honcnlty.  ,hlt  ls 
unimportant,  as  long  as  th.  isolated  control  surface  is  not  used 
for  accurate  observations  of  earth  tides.  This  proble.  can  l. 
elir.in.ted  by  oaklne  the  isolated  control  surface  completely 
closed  except  for  thensally  insulated  electrical  feed  tin.- 


Th»  open  loop  g«in  of  the  temperature  control  system 
extrea^  high  below  10**  Ht  end  fell,  70  db/dec  above  thi,  fre¬ 
quency.  This  active  device  provide,  good  regulation  where  the 
passive  Insulation  is  ineffective.  Conversely,  \he  insulation 
attenuates  teapereture  fluctuation,  at  higher  frequencies,  where 
the  low  loop  gain  reduces  the  effectiveness  of  the  teapereture 
controller.  A  phase  lead  coapensator  is  use<i  to  stabil¬ 

ize  the  system  if  needed.  A  series  of  tests  of  the  temperature  con¬ 
troller  with  varying  thermal  time  constants  between  the  cont^T  heater 
and  the  outer  aluminum  can  show  that  from  approximately  5  cph  to  the 
seisaic  band,  no  appreciable  noise  comet  from  the  controller. 

Oie  possible  source  of  thermally  Induced  noise  within 
the  mechanical  instrument  Itself  Is  ohmic  heating  in  the  center 
transducer  plate  and  auxiliary  damping  paddles  due  to  the  eddy 
currents  which  provide  the  damping.  If  unwanted  modes  are 
variably  excited,  this  heating  could  conceivably  caus*  fluctua¬ 
tions  in  the  fibrr  temperature.  A  worst  case  analysis  of  this 
effect  proceeds  as  follows.  Suppose  that  one  of  the  unwanted  modes 
of  vibration  in  the  vertical  accelerometer  is  continually  excited 
»y  s*te  noise  to  an  amplitude  corresponding  to  full  scale  open 
loop  output.  Then  the  energy  in  the  mode  is  approximately 


where  m  is  the  accelerometer  mass,  «*2.f  i,  the  angular  frequen- 
V  of  the  mode,  and  6fg  i,  the  full  scale  displacement  of  the  mass. 
If  0  is  the  quality  factor  for  the  mode,  the  amount  of 
energy  dissipated  per  cycle  is  by  definition  of  Q,  Thus,  the 
power  lost  to  heating  can  be  written  as 


P  =  2lLf  e  £±»  «  2 

p  Q  20*“  6fs 


(2) 


with  the  help  of  equation  (1). 

frequencies  end  Q's  for  typical  unwanted  modes  in  the 
vertical  accelerometer  have  been  estimated  as  follows: 

f  •  40Hz,  Q  •  30  (3) 

The  parameters  m  and  are  roughly  lCg  and  5xl<T*cm, 

respectively.  Hence,  the  rate  of  heat  production  in  the  dashing 
plates  due  to  continuous  excitation  of  an  unwanted  mode  is 

.5  ,tj 

P  7  6x10  erg/sec  »  1.5x10  cal/sec  (4) 

Since  the  power  (equation  2)  depends  on  the  cube  of  the  modal 
frequency  u,  the  heating  due  to  continual  excitation  of  the 
fundamental  mode  would  be  considerably  less  than  this. 

Now  the  heat  capacity  of  the  entire  mass  assembly  is 
roughly  that  of  the  copper  mass  itself,  or  about  Ch  *  2cal/°C. 
If  all  the  energy  dissipated  in  damping  is  assumed  to  heat  this 
assembly,  then  the  rate  of  change  of  its  temperature  is 

=  P/C.  *  l0'12°C/sec  (5) 

at  h 

This  rate  is  just  large  enough  to  cause  a  barely  perceptible 
temperature  change  of  10  °C  in  one  day  of  continuous  excitation 
of  an  unwanted  mode.  Moreover,  heating  due  to  damping  of  the 
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fundaoental  mode  is  «uch  slower.  Obviously,  ths  offset  of  heating 
in  the  daaping  assemblies  is  negligible  ss  far  as  the  accelerometer 
output  is  concerned. 
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filtering  and  Recording  of  Pat* 

The  quartz  fiber  accelerometers  in  the  open  loop  systea  have 
broadband  response  from  DC  to  the  resonance  frequoncy  of  about  1  Hz. 
with  this  bandwidth,  simultaneous  investigations  of  earth  tides,  frs*e 
oscillations  of  the  earth,  and  long  period  teleseisms  can  be  carried  out. 
Tidt.1  data  is  sore  easily  analyzed  when  earthquakes  are  removed  froo  it. 
Teleaeismic  records  can  be  obtained  at  a  higher  gain  when  the  tidal  signals 
are  removed.  Earth  normal  mode  studies 6are  conveniently  studied  through  a 
band  pass  gainful  filter.  The  raw  output  of  the  accelerometer  Is  thus 
ro°dlfi*d  by  thr®«  ««ive  analog  filters  called  prefilters,  which  respond 
preferentially  to  tides,  free  oecillatlone,  and  teleseisms.  The  resultant 
total  response  functions  are  plotted  along  the  accelerometer  response  for  the 
cpen  loop  systea  in  figure  6. 

Each  of  the  three  prefilters  is  an  electronic  circuit  which  uses 
fedback  operational  amplifiers  to  set  Its  response  and  the  full  scale 
output  of  ♦  10V.  The  tidal  prefilter  simply  attenuates  high  frequency 
signals,  while  the  normal  rode  and  se.snic  prefilters  both  have  gainful 
response  at  certain  frequencies. 

The  tidal  prefilter  is  a  double-pole  lowpass  filter  with  a  cutoff 
frequency  of  5.75  cph  (slow)  or  57.5  cph  (fast).  Tor  tidal  investigations, 
the  lower  cutoff  frequency  is  used.  Since  the  highest  tidal  frequency  of 
usual  Interest  is  that  of  the  triurnal  tide,  or  0.175  cph,  the  slow  cutoff 
frequency  gives  the  tidal  prefiltcr  nearly  unity  gain  and  zero  phase  shift 
over  the  entire  tidal  band. 

The  response  near  DC  is  used  for  calibration  of  the  instrument. 

The  tidal  prefilter  is  used  during  calibration  to  observe  DC  levels  while 
averaging  out  fluctuations  due  to  aicroseismic  ground  noise.  To  avoid  long 
waiting  periods,  the  time  constant  is  set  to  fast  for  such  work. 


Normal  mode  oscillations  of  the  aarth  arc  best  investigated  using  a  gainful 
bandpass  filter.  Its  passband  is  variable,  with  a  lower  knee  frequency  of 
1  eph  and  an  upper  knee  frequency  of  30,  60,  120  or  180  cph.  The  gain  in 
the  passband  is  40db  above  the  DC  level.  On  the  high  side  of  the  passband, 
the  response  falls  40  db/dec.  On  the  low  side,  it  falls  40db/dec  to  0.1  cph 
and  is  constant  at  lower  frequencies.  An  input  attenuator  sets  the  DC  level. 
Vhen  the  attenuation  is  zero,  the  gain  at  DC  is  unity. 

The  normal  mode  prefilter  is  usually  run 
with  no  input  attenuation  and  the  upper  knee  frequency  set  at 
30  cph.  This  configuration  provides  sufficient  gain  and  bandpass 
for  studies  of  the  lower  frequency  normal  inodes.  Increasing  the 
upper  knee  frequency  includes  more  of  the  normal  node  band  in  the 
gainful  region.  At  the  highest  setting,  long  period  seismic  signals 
are  included  in  the  bandpass.  This  setting  allows  all  seismic 
signals  with  periods  greater  than  20  sec  to  be  observed  with  near!; 
the  same  gaim.  Unfortunately,  larg*  Rayleigh  waves  and  ground  noise 
tend  to  saturate  the  prefilter  electronics  when  its  bandwidth  is 
this  large.  This  is  the  reason  the  input  attenuator  Is  provided. 

It  allows  a  compromise  between  gain  and  bandwidth  to  be  made.  The 
nature  of  this  compromise  depends  on  the  size  of  events  to  be 
observed  and  the  level  of  ground  noise  at  the  instilment  site. 

The  normal  mode  prefilter  attenuates  the  t  .dal  signals  which  allows 
more  gain  to  be  used  in  the  bandpass.  Peak-to-peak  tidal  accelera- 
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tions  can  be  as  high  as  3xl0~7g,  while  normal  modes  and  seismic 
waves  typically  produce  accelerations  of  order  10  *g.  Tidal  signals 
at  the  lock-in  amplifier  output  are  thus  near  IV  peak-to-peak.  If 
amplified  significantly,  these  signals  saturate  the  prefilter. 

Phase  shifts  due  to  the  changing  gain  of  the  prefilter  near  tidal 
frequencies  distort  the  tidal  signals.  This  precludes 

use  of  the  normal  mode  prefilter  output  for  tidal  studies  or  tidal 
calibration. 

The  seismic  prefilter  is  a  peaked  bandpass  filter.  Its  peak 

is  fixed  at  100  second  periods,  while  its  half  power  points  are  at 
60  and  200  second  periods.  At  high  frequencies,  its  response  falls 
as  w'7  ,  or  140db/dec,  while  at  low  frequencies  its  response  dies 
as  u2  ,  or  40db/dec.  The  peak  gain  is  set  by  the  input  attenuator, 
which  has  six  positions,  and  can  be  varied  from  60  to  5830.  Normally, 
the  attenuator  is  set  for  a  peak  gain  of  621. 

The  high  peak  gain  and  the  rather  narrow  bandwidth  of  the 
seismic  prefilter  are  quite  useful  in  observing  long  period  tele- 

-10 

seisms.  They  permit  measurement  of  seismic  signals  as  low  as  10  g 
in  amplitude.  The  sharp  rolloff  of  the  response  at  high  frequencies 
reduces  interference  due  to  microseismic  signals,  while  the 
rolloff  at  low  frequencies  eliminates  the  tides.  The  seismic  pre¬ 
filter  is  used  primarily  for  detection  of  distant  teleseismic  events. 

The  accelerometers  have  flat  responses  to  acceleration 
over  the  bandpasses  of  the  prefilters.  Thus,  the  shape  of  the 
composite  response  at  the  output  of  each  p>refilter  is  determined 
by  the  response  of  the  prefilter  itself.  The 
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gain  level  ia  determined  by  the  DC  gain  of  the  accelerometer.  The 
exact  form  of  the  response  to  acceleration  at  each  prefilter  output 
is  calculated  by  multiplying  the  prefilter  response  function  by  the 
appropriate  accelerometer  response. 

The  output  of  each  prefilter  is  recorded  digitally  on  magnetic 
tape  with  a  resolution  of  +16384  digitizer  counts  per  +10  volts. 

In  the  cases  of  the  normal  mole  and  seismic  prefilters,  the  outputs 

with  +10V  full  scale,  are  used  for  digitizing.  The  typical  sampling 
interval  is  1  second,  giving  a  digitizer  Nyquist  frequency  of  1/2  Hz. 
Part  of  the  function  of  the  prefilters  is  to  provide  large  attenuation 
above  the  Nyquist  frequency  to  eliminate  aliasing  of  energy  in  the 
digitizing  process.  As  Figure  6  shows,  the  response  of  each  prefilter 
at  the  Nyquist  frequency  is  down  at  least  80  db  from  its  peak  response. 
Such  attenuation  is  sufficient  to  eliminate  aliasing  for  all  practical 
purposes. 
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Noise 


Measurements  cf  the  fuU  electronic  system  noise  „ill  be  presented 
below  end  compared  with  calculated  Brownian  motion  noise  and  measured 
ambient  site  noise.  This  will  be  done  for  the  normal  mode  and  seismic 
prefilter  channels.  The  calculated  Brownian  noise  of  the  mass-torsion 
fiber  system,  with  appropriate  measured  gain  factors  and  measured  pre¬ 
filter  response  functions  is  contained  in  Table  I  as  an  equivalent  output 
voltage.  The  electronics  is  assumed  to  produce  no  noise  in  these  calculations. 

The  full  electronic  system  noise  is  measured  with  the  instrument 
in  thermal  control  and  with  the  synchronous  detector  drive  voltage  set 
to  zero.  The  measured  values  are  contained  in  Table  II  for  the  normal  mode 
prefilter  and  the  seismic  prefilter  out  :s.  Comparison  of  the  values  in 
Tables  I  and  II  show  that  the  noise  in  the  normal  mode  prefilter  output  from 
both  sources  is  about  equal  and  the  electronic  noise  is  lower  than  the 
Brownian  noise  in  the  seismic  prefilter  output. 

A  typical  ambient  output  in  the  seismic  filter  is  shown  in  rig.  7 
for  both  day  and  night.  It  should  be  noted  that  there  is  little  difference 
in  these  levels.  The  peak  to  peak  output  voltage  is  about  500  mv  which  is 
considerably  larger  than  the  combined  Brownian  and  electronic  noise  levels. 

The  microseWarejisible  in  both  traces  of  Fig.  7  as  high  frequency  hash. 

The  largest/wiggles  are  at  about  100  seconds ,  the  period  at  which  the 
seismic  prefilter  response  is  peaked.  Recent  studies8l9'11,have  shown  that  the 
effects  of  propagation  of  pressure  cells  in  the  atmosphere  become  important 
at  periods  greater  than  60  seconds.  These  wiggles  are  probably  due  to 
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local  ground  motion  induced  by  these  pressure  cells.  The  presence  of 
this  noise  just  below  the  seismic  band  (20  -  50  second  periods)  makes  the 
use  of  digital  filtering 5an  important  technique  to  improve  the  signal  to 
noise  ratio  for  teleseismic  events.  The  ambient  noise  in  the  normal  mode 
prefilter  with  the  upper  knee  frequency  of  30  cph  and  no  input  attenuator  was 
between  50  and  100  mv  peak  to  peak.  A  comparison  of  Tables  I  and  II  show 
that  the  ambient  noise  is  larger  than  the  combined  electronic  and  Brownian 
noise. 

The  improvement  in  noise  performance  of  the  present  instrument 
compared  to  the  first  model  of  the  vertical  accelerometer  is  between  6  to 
10  db  in  the  normal  mode  (1-30  cph)  prefilter  band.  The  decrease  in  noise 
lor  the  seismic  prelilter  band  is  about  a  factor  of  three  in  amplitude. 
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Data 


The  modified  instrument  described  above  has  been  in  operation 
since  February  1971.  During  this  time  a  continuous  digitally  recorded 
run  was  made  at  Camp  Elliott,  San  Diego  from  February  to  July  inclusive 
to  obtain  data  from  a  sizable  number  of  earthquakes.  This  data  was  taken 
with  the  prefilters  described  in  the  section  Filtering  and  Recording  of  Data. 

A  complete  catalogue  of  the  events  detected  is  in  progress  by  T.  Whorf  at  IGPP. 

Selected  examples  of  earthquakes  as  seen  through  the  normal  node 
prefilter  and  seismic  prefilter  are  given  below.  Fig.  8  shows  the  normal 
mode  prefilter  record  of  an  earthquake  =  6.3)  in  the  Southern  Shetland 
Islands  (63.5°S,  61.2°W;  on  February  8,  1971,  origin  time  21:04  GMT  at  a 
depth  of  33  km.  The  ambient  noise  through  the  prefilter  before  the  earth¬ 
quake  is  also  seen  in  Fig.  8.  The  power  spectrum  of  this  earthquake  from 
15  hours  of  data  is  plotted  in  Fig.  9  where  zero  db  corresponds  to  1>;10~*  *g^/chp. 
The  mean  of  the  ambient  noise  power  spectrum  is  also  plotted  on  Fig.  9. 

The  ambient  noise  and  its  mean  are  shown  displaced  downward  from  the  r  true 
values  by  10  db  on  Fig  9  so  as  not  to  be  confused  with  the  earthquake  power 
spectrum.  The  ambient  noise  power  spectrum  is  derived  from  the  same  length 
of  record  used  in  the  earthquake  power  spectrum  and  is  processed  by  the  same 

program.  In  this  case,  data  from  the  15  hours  proceeding  the  earthquake 
was  used. 

The  time  record  of  this  earthquake  (Fig.  8)  shows  several  passes 
of  the  Rayleigh  waves  and  commensurately  the  power  spectrum  shows  clearly 
defined  normal  mode  excitation  with  a  20  db  or  larger  signal  to  noise  ratio 
above  6  cph.  Unfortunately  this  fine  record  was  interrupted  bv  the 
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San  F ernando  earthquake  (H^  =  6.2)  on  February  9,  1971,  origin  time  14:01  GMT. 
Tigs.  10  and  11  show  the  normal  mode  prefilter  record  of  this  event.  Since 
the  accelerometer  was  essentially  at  the  epicenter .instrumental  saturation 
occurred  at  the  primary  event  and  for  the  larger  aftershocks.  There  are 
two  main  features  to  be  noted  in  rigs.  10  and  11.  Being  at  the  epicenter 
means  that  both  Rayleigh  wave  packets  arrive  simultaneously  and  the  long 
wave  length,  high  velocity  Rayleign  waves11  can  be  seen  with  minimal  inter¬ 
ference  from  the  Rayleigh  wave  packets.  The  aftershocks  seen  in  rigs.  10 
and  11  should,  by  their  form  and  magnitude,  yield  information  on  the  moment 
of  the  event  and  its  source  structure.  An  investigation  of  the  aftershocks 
will  be  made  for  this  purpose.  The  tidal  prefilter  shows  less  saturation 
than  the  earth  mode  prefilter  because  of  its  lower  gain.  The  tidal  pre¬ 
filter  can  provide  better  information  on  the  aftershocks  for  thin  reason. 

No  net  DC  displacement  in  the  tidal  channel  was  observed  larger  than  3xl0~^®g 
or  an  equivalent  vertical  motion  of  1  ran  relative  to  the  center  of  mass  of 
the  earth.  This  estimate  is  an  uppir  bound  only  owing  to  difficulties  in 
matching  the  endpoints  of  the  tidal  signal  because  of  saturation  problems. 

Fig.  12  shows  an  earthquake  (Mfa  =  6.0)  in  Adak  (51.4®N,  177. ?°K) 

•n  May  1971  origin  time  06:06  GMT  at  a  depth  of  43  km.  The  pewer  spectrum 
of  this  earthquake  from  18  hours  of  data  is  plotted  in  Tig  14.  The  ambient 
noise  and  its  mean  are  plotted  as  in  Tig.  9.  The  ambie  it  noise  is  obtained 
from  the  previous  18  hours  of  normal  mode  prefilter  lata.  The  power  spectrum 
of  -his  earthquake  is  quite  different  from  that  of  the  Shetland  Island  event. 

he  underlying  reasons  for  the  general  features  of  these  power  spec-ra  is  &,!il 
not  clearly  understood  although  computation  of  excitation  spectra  from  earth 
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quake  nodels  Is  now  possible.  A  spectrum  timilar  to  those  shown  In  Tigs.  9 
and  11  is  obtained  somewhat  more  than  once  a  month  on  an  average. 

Seiseic  prefilter  data  is  contained  in  Tigs.  13  and  l*i  where  a 
typical  series  of  small  teleseismic  events  taken  over  several  days  are 
shown.  Table  II  contains  identifications  and  pertinent  data  on  the  events. 

il 

Kvent  17  has  been  previously  published  and  is  omitted.  The  salient  feature 
of  this  data  is  the  low  ambient  noise  level  through  the  seismic  prefilter 
and  the  consequent  low  magnitude  threshold  for  visual  detection  of  events. 
The  visual  detection  threshold  can  be  greatly  improved  by  digital  high  pass 

4  t  s 

or  band  pass  filtering.  The  seismic  prefilter  data  shown  is  typical  in 
its  noise  performance  for  the  entire  body  of  data  taken  from  February  to 
July. 

A  paper  on  the  horizontal  analogue  of  this  vertical  accelerometer 
design  with  preliminary  data  is  in  progress. 
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TABLE  I 
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TABLE  II 


Normal  Node  Prefilter 
with  no  input  attenuation 

Upper  Knee  Frequency 

30  cph 


Seismic  Prefilter 
with  peak  gain  at  621 
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FIGURE  CAPTIONS 


Figure 

Figure 
Figure 
Figure  4 

Figure  5 

Figure  6 

figure  7 


Assembly  of  the  mess-torsion  fiber  system  with  the  drive  capacitor 
plates  in  place. 

Exploded  view  of  the  fiber  center  clamp  assembly  with  mass. 

Assembled  fiber  center  clamp  assembly  with  mass. 

Assembled  mass-torsion  fiber  system  with  assembly  of  mode  damping 
magnets.  The  pitted  surface  of  the  aluminum  frame  comes  from 
chemical  cleaning. 

Thermal  insulation  system  for  the  vertical  accelerometer.  Note 
the  temperature  controller  heater  windings  on  control  surface. 

Composite  responses  to  acceleration  of  open  loop  instrument  and 
three  prefilters,  with  response  at  position  sensor  output  shown 
for  comparison.  All  response  magnitudes  are  plotted  relative  to 
open  loop  DC  gain  of  accelerometer  which  is  6.24xl06  volts/g. 

Daytime  and  nighttime  vertical  noise  in  seismic  band  at  Elliott 
Station.  Records  are  computer  plots  of  digitized  output  of 
seismic  prefilter  for  cpen  loop  vertical  accelerometer.  Input 
attenuator  of  prefilter  is  set  at  a  gain  of  621.  Scales  apply 
to  both  traces.  Calibration  in  terms  of  ground  acceleration  can 
be  made  using  Figure  7  with  DC  gain  taken  as  6.24xl06  volts/g. 
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figure  8 


figure 


figure 


figure 


Coaputer  plot  of  normal  node  prefilter  data  for  Southern 
.Shetland  Island  earthquake  (M^  *  6.3).  Vertical  scale  corresponds 
only  to  signals  in  1-30  cph  band-pass,  not  to  tides.  The 
epicenter  was  <63.5«S,  61.2®W),  origin  tine  21:04  GMT,  Tebruary  8, 
1971  at  a  depth  of  33  km. 

9  Spectra  of  the  Southern  Shetland  Island  earthquake  (Tig.  8) 

from  2  to  30  cph.  Absolute  spectral  density  of  0  db  point  is 
-24  2 

1x10  g  /cph.  Earth  normal  mode  spectrum  is  derived  from  15 
hours  of  data  shown  in  fig.  8.  Ambient  spectrum  is  derived  from 
15  hours  of  data  just  previous  to  the  earthquake.  The  running 
mean  of  ambient  spectrum  is  average  over  1  cph.  The  ambient 
spectrum  and  running  mean  are  displaced  downward  by  10  db  In 
Fig.  9.  Approximate  level  of  Brownian  noise  at  input  of  accelero¬ 
meter  it;  indicated  by  heavy  line.  Selected  modal  lines  arc 
identified  from  theoretical  calculations  of  eigenfrequencies. 

10  Computer  plot  of  normal  mode  prefilter  data  for  the  San  Ternando 

earthquake  s  6.2).  Vertical  scale  corresponds  only  to  signals 

in  1-30  ph  bandpass,  not  to  tides.  The  origin  time  was  14:01  u.TT, 
February  9,  1971,  and  the  epicenter  was  (34.4®N,  118. 4°W)  at  a 
depth  of  13  km. 

11  Continuation  of  figure  10. 
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Figure  12 


Tiguro  13 


figure  14 


Figure  15 


Computer  plot  of  normal  node  prefilter  date  for  the  Adak 
earthquake  (M^  *  6.0).  Vertical  scale  corresponds  only  to 
signals  in  1-30  cph  bandpass,  not  to  tides.  The  epicenter  was 
(51.4°K,  177. 2°W),  origin  time  06:08,  May  2,  1971  at  a  depth 
of  43  km. 

Spectra  of  the  Adak  earthquake  (rig.  12)  from  2  to  30  cph. 

Earth  normal  mode  spectrum  is  derived  from  18  hours  of  data  shown 
in  Tig.  12.  the  ambient  spectrum  is  derived  from  18  hoars  of 
data  just  previous  to  the  earthquake.  See  caption  for  rig.  9 
for  explanation  of  ambient  and  mean  spectra. 

Computer  plot  of  seismic  prefilter  data  for  a  series  of  3itall 
teleseisric  events.  (See  Table  II  for  pertinent  data.)  Event 
17  has  been  published  (see  text). 

See  caption  figure  14 
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ABSTRACT 


A  wide  band  horizontal  accelerometer  using  a  quartz  torsional 
pivot  and  capacitative  position  sensing  with  phase-sensitive  detection  is 
described.  Devices  for  magnetic  damping  of  unwanted  modes  of  vibration 
and  internal  stabilization  of  temperature  and  pressure  are  discussed  in 
detail.  The  accelerometer  is  small  enough  in  size  and  weight  to  allow 
deployment  in  bore  holes.  Data  from  the  earth  normal  mode  band  and  the 
seismic  band  are  shown.  Although  horizontal  noise  due  to  surface  tilting 
is  high  for  periods  greater  than  50  seconds,  noise  in  the  seismic  band  can 
be  reduced  by  as  much  as  an  order  of  magnitude  by  digital  high  pass  filter¬ 
ing.  The  resulting  horizontal  signal-to-noise  ratio  for  teleseismic  events 
is  comparable  to  that  for  the  vertical  signals. 


1 


6  7 


1.  The  Pendulum 


Due  to  the  presence  of  gravity,  sensitive  vertical  accelero¬ 
meters  are  m  general  harder  to  make  than  their  horizontal  counter¬ 
parts.  In  vertical  instruments,  some  spring  device  must  precisely 
cancel  the  acceleration  of  gravity.  Since  the  residual  after 
cancellation  must  be  below  the  desired  level  of  measurement,  strin¬ 
gent  stability  requirements  are  imposed  on  the  spring  device.  Any 
changes  in  the  spring  force  which  opposes  gravity  appear  as  spurious  out¬ 
put.  Fortunately,  this  problem  does  not  occur  in  horizontal  measure¬ 
ments.  Nearly  all  horizontal  accelerometers  are  variants  of  two 
basic  devices:  the  pendulum  and  the  "swinging  gate."  Stability 
problems  do  occur  in  these  systems,  as  the  pivots  or  supports  must 
work  smoothly  and  consistently  over  long  time  periods.  However,  the 
problem  here  is  not  one  of  stability  under  stress. 

Many  different  suspension  systems  have  been  used  in  hori¬ 
zontal  accelerometers.  Among  them  are  knife  edges,  cross-spring 
pivots,  and  magnetic  supports. H ,16,17  At  th(?  infinitesimal  levels 
of  motion  and  force  involved  in  geophysical  work,  most  of  these 
devices  have  instabilities.  They  often  pop  or  creak,  and  pivots 
made  of  metal  tend  to  drift.  The  present  design  takes  advantage 
of  the  excellent  stability  of  a  quartz  fiber  suspension  developed 
for  use  in  a  vertical  accelerometer.6  The  fiber  suspension 
is  made  from  a  fused  quartz  rod  0.19  inches  in  diameter  with  ground 
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clamping  flats  at  Its  center  and  ends.  Two  active  sections,  more 
than  ten  times  smaller  in  diameter  and  each  about  2.5  cm  long,  are 
drawn  between  the  center  and  end  sections.  This  design  eliminates 
instabilities  in  the  clamps,  a  common  source  of  trouble  -in  fiber 
suspensions,  by  providing  large  flats  for  positive  contact. 

At  the  fiber's  center,  a  pendulum  assembly  is  suspended. 
The  pendulum  is  preferred  to  the  "swinging  gate"  design,  as  the 
latter  has  a  fundamental  asymmetry  with  respect  to  gravity.  The 
mass  of  the  pendulum  is  a  flat  plate  which  aids  in  the  magnetic 
damping  of  unwanted  modes  of  vibration.  It  is  sufficiently  large, 
and  the  fiber  is  sufficiently  weak,  so  that  gravitational  restoring 
forces  predominate.  Thus  the  fiber  serves  only  as  a  stable 
loss-free  pivot,  and  the  free  period  of  the  device  is  its  pendulum, 
period.  Motion  of  the  pendulum  is  detected  by  means  of  a 
three-plate  capacitative  position  sensor,  whose  center  plate  is 
attached  to  the  pendulum  assembly.  The  pendulum  length  is  75  cm. 
During  construction,  the  fiber  is  rotated  so  that  it  is  free  of 
torsion  when  the  pendulum  hangs  vertically.  However,  the  fiber  is 
kept  under  tension  by  means  of  an  elinvar  parallel-motion  spring. 

The  position  sensor  and  the  mechanisms  for  tensioning  and  rotating 
the  fiber  are  the  same  as  those  in  the  vertical  accelerometcr.6 
figure  1  shows  the  fiber  suspension  arid  pendulum  assembly. 

To  avoid  direct  magnetic  forces  on  the  moving  parts,  all 
materials  used  in  the  pendulum  assembly  are  free  from  magnetic 
contamination.  Where  strength  is  necessary,  beryllium  copper  is 
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FIGURE  1 


Photograph  of  fiber  suspension  and  pendulum  assembly  for  horizontal 
accelerometer.  Fiber  center  section  and  active  sections  are  visible 
in  center  of  frame.  Large  detector  plate  serves  as  center  element 
of  three  plate  capacitative  position  sensor.  Small  plates  attached 
a  fiber  clamp  are  auxiliary  damping  paddles,  which  serve  to  reduce 
excitation  of  unwanted  vibrations.  One  of  lower  auxiliary  damping 
magnets  is  visible  below  paddle.  Plate  at  right,  in  conjunction  with 
elinvar  parallel  motion  spring,  maintains  tension  in  fiber,  which 
is  set  with  large  bolt  at  right.  Slip  rings  under  fiber  end  clamps 
allow  rotation  of  fiber  during  assembly.  Frame  housing  assembly 
mates  with  similar  frame  supporting  outer  plates  of  position  sensor 
(of.  Figure  2).  Tripod  is  not  part  of  instrument. 


FIGURE  2 


Position  Sensor  Assembly.  Note  the  damping  magnets  on  the  capacitor 
plates  and  plink  stoppers  in  the  arc  of  the  magnets.  Plates  shown  are 
from  prototype  vertical  accelerometer.  Those  actually  used  for 
horizontal  are  twice  as  wide. 
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used.  Otherwise,  nearly  pure  (1100  series)  aluminum  or  OFHC 
copper  is  used.  Even  the  tiny  0-80  screws  in  the  assembly  are 
non-magnetic,  being  fashioned  from  beryllium  copper. 

Figure  3  ,  shows  an  exploded  drawing  of  the  pendulum 

assembly.  Beryllium  copper  clamps  B6  and  B9  mate  with  the  ground 
flats  on  fiber  center  section  B7.  Pendulum  mass  B12,  cf  OFHC 
copper,  is  supported  by  beryllium  copper  rod  Bll,  which  is  con¬ 
tracted  in  the  drawing.  Aluminum  detector  plate  Bl,  part  of  the 
position  sensor,  is  attached  to  clamp  B6  by  means  of  quartz  insu¬ 
lator  B5  and  beryllium  copper  clamping  bars  B2-B4.  The  purpbse 
of  aluminum  auxiliary  damping  paddles  B8  and  BIO  and  pendulum 

damping  magnets  B13  and  B14  is  to  aid  in  reducing  excitation  of 
unwanted  modes. 

In  the  fundamental  mode  of  vibration,  the  pendulum 
executes  small  rotational  oscillations  about  the  fiber  axis.  For 
a  small  rotation  Q  away  from  the  vertical,  the  restoring  torque 
due  to  gravity  has  the  form 

tg  -  -mLmg  sine  =  -mLmge  (lll) 

where  m  is  the  mass  of  the  entire  pendulum  assembly,  including 
the  support  and  clamps,  and  Lm  is  the  distance  of  its  center  of 
mass  from  the  fiber  axis.  This  restoring  torque  corresponds  to 
an  effective  gravitational  torsion  constant  mL^g.  If  the  fiber 
has  a  torsion  constant  K  and  the  pendulum  assembly  has  a  moment 
of  inertia  I  about  the  fiber  axis,  the  equation  of  undamped 
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FIGURE  3 


Exploded  diagram  of  pendulum  assembly  for  horizontal  accelerometer. 
Fiber  center  section  B7  is  held  by  fiber  clamps  B6  and  B9  ,  and 
pendulum  mass  B12  is  supported  by  rod  Bll.  Detector  plate  B1  is 
is  electrically  insulated  from  rest  of  assembly  by  quartz  insulator 
B5.  Permanent  magnets  B13  and  B14  damp  unwanted  vibrations  by 
inducing  eddy  currents  in  pendulum  mass  B12.  They  are  supported 
by  sleeve  (not  shown)  surrounding  rod  Bll.  Plates  B8  and  BIO  are 

auxiliary  damping  paddles.  Auxiliary  damping  magnets  are  not 
shown. 


motion  is 


IB  =  t  -  K6 

£ 

or  18  =  -<mLmg  ♦  K)0  (1.2) 

..  2 

or  0  t  co  6  =o 

O 

where  +  ^ 

o  — - 

Equation  (1.1)  has  been  used  for  the  gravitational  restoring 
torque.  Thus  for  small  amplitude  0  the  fundamental  mode  is  an 
harmonic  oscillation  with  a  free  period 

T°  2  —  '  f/CET  a  3) 

%  /  mL  g  +  K  (1.3) 

u  m 

The  mechanical  parameters  of  the  instruments  have  the 
following  approximate  values 

m  =  35  g 

Lm  =  15  cm  (1.4) 

I  =  1.36  x  10^g-cm^ 

Ihsse  parameters,  the  effective  torsion  constant  of 

gravity  is  mL  g  =  5.1  x  105  dyne-cir 

rad 
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In  the  vertical  instrument , ’the  fiber,  which  is  torsionally  stressed 
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to  support  the  weight  of  the  mass,  has  a  torsion  constant  of 
7.5x10^  - .  The  torsion  constant  of  the  unstressed  fiber 
in  the  horizontal  must  be  even  less.  Thus,  th»*  'iber  in  the 
horizontal  supplies  at  most  1.5%  of  the  restoring  torsion.  As 
stated  above,  the  horizontal  accelerometer  acts  like  a  simple 
pendulum.  Its  free  period  is  given  by  equation  (1.3),  with  K 
being  ignored  as  a  small  quantity.  Substitution  of  the  parameters 
(1.4)  yields  T^  =  1.0  sec.  The  measured  natural  period  is 


Tq  =  0.910  sec 


(1.5) 
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2.  Unwanted  Modes  of  Vibration 


The  ideal  mechanical  system  for  an  accelerometer  is  one 
with  strictly  a  single  degree  of  freedom.  Unfortunately,  this 
ideal  cannot  be  physically  realized,  as  every  real  system  is 
three-dimensional  and  so  has  more  than  one  mode  of  vibration. 
Residual  non-linearities  transfer  energy  from  the  unwanted  modes 
to  the  fundamental  mode.  This  coupling  can  be  a  major  source  of 
noise  when  the  unwanted  modes  are  highly  excited.6 

The  fiber  suspension  has  two  advantages  over  other  elastic 
supports  in  this  respect.  It  has  a  high  degree  of  symmetry,  and 
it  has  an  extra  variable,  the  fiber  tension,  with  which  to  control 
the  unwanted  modes.  In  the  limit  of  infinite  tension,  the  un¬ 
wanted  modes  are  absent  and  the  fiber  can  only  rotate.  As  shown 
by  equations  (1.2)  ,  the  undamped  frequency  of  the  fundamental 
mode  is  independent  of  the  fiber  tension.  All  of  the  unwanted 
modes,  however,  have  frequencies  which  increase  as  the  fiber  ten¬ 
sion  is  raised.  The  fiber  tension  is  adjusted  so  that  the 
frequencies  of  the  unwanted  modes  are  separated  from  that  of  the 
fundamental  and  placed  in  regions  of  low  seismic  noise.  This 
protects  the  fundamental  mode  from  non-linear  coupling.  Additional 
protection  is  provided  by  magnetic  damping  assemblies  which  affect 
the  unwanted  modes  preferentially. 

Besides  the  fundamental  mode,  the  fiber-pendulum  system 
has  four  significant  normal  modes  of  oscillation,  two  translational 
modes  and  two  twisting  modes.  The  translational  modes  involve 
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horizontal  and  vertical  motion  of  the  fiber  center  perpendicular 
to  the  fiber  axis,  while  the  twisting  modes  involve  rotation  of 
the  pendulum  assembly.  In  one  twisting  mode  the  assembly  rotates 
about  its  vertical  axis  of  symmetry.  In  the  other,  it  rotates 
about  a  horizontal  axis  perpendicular  to  the  fiber.  Together  with 
the  fundamental,  the  four  modes  comprise  the  three  principal 
rotations  and  two  of  the  principal  translations  of  the  system. 
Translation  along  the  fiber  axis  is  precluded  by  the  supporting 
structure . 

Crude  estimates  of  the  frequencies  of  the  unwanted  modes 
can  be  made  from  simple  physical  reasoning.  For  this  purpose  the 
pendulum  assembly  is  characterized  by  its  total  mass  m  and  its  moments 
of  inertia  1^  and  I^  about  the  two  twist  axes.  It  rests  on  the 
fiber  center  section  of  length  £q,  which  is  supported  by  the  two 
active  sections  of  length  i  under  tension  T.  Figure  4  illus¬ 
trates  the  type  of  motion  involved  in  each  mode. 

The  two  translational  modes  are  functionally  identical. 

As  the  center  section  of  the  fiber  translates  a  small  amount  x, 
the  active  fiber  sections  make  angles  0  =  x/£  with  their  relaxed 
directions.  The  active  sections  are  assumed  to  remain  straight 
and  unchanged  in  length,  and  their  stiffness  is  neglected.  Under 
these  assumptions,  the  restoring  force  is  provided  by  the  component 
of  the  tension  in  the  direction  of  the  translation.  Since  there 
are  two  active  fiber  sections,  the  restoring  force  is 
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FIGURE  4 


Schematic  illustration  of  unwanted  modes  of  vibration  of  pendulum 
assembly  in  horizontal  accelerometer.  Active  fiber  sections  of 
length  £  are  assumed  neither  to  stretch  nor  to  flex.  Fiber 
center  section  has  length  t0.  Entire  pendulum  assembly  has  mass 
m  and  moments  of  inertia  1^  ^  about  two  twist  axes 

(a)  Vertical  translation  mode. 

(b)  Horizontal  translation  mode. 

(c)  First  twisting  mode.  Rotation  is  about  vertical 
axis  of  symmetry  of  pendulum  assembly. 

<d>  Second  twisting  mode .  Rotation  is  about  horizon- 
tal  axis  perpendicular  to  fiber  (out  of  paper  in 
drawing). 
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F  =  -2T  sin0  =  -2T0  =  -  27x  (  2.1) 

1 

The  equation  of  notion  for  the  translations  is 


M 

ox  = 

F  =  -2T 

X 

I 

or 

X  + 

2 

«t  x  =  0 

where 

2 

"t  8 

2T 

mi 

Thus  the  translational  modes  have  a 

common  frequency 

f*  = 

Ut  *  JL__ 

/TT 

t 

2w  2»r 

/  ml 

As  promised,  this  frequency  increases  with  the  tension  in  the 
fiber,  although  only  as  its  square  root. 

The  two  twisting  modes  are  nearly  alike,  tr.^ir  only 
difference  being  in  the  moments  of  inertia.  As  the 
pendulum  assembly  twists  through  a  small  angle  Q  ,  the  two  active 
fiber  sections  make  angles  O'  with  their  original  directions. 
Again  their  lengths  are  assumed  to  remain  unchanged  and  their 
stiffness  and  bending  is  ignored.  From  the  geometry  shown  in 

Figure  4  ,  the  following  relation  between  the  two  small  angles 
0  and  O'  is  obtained: 


16 


9t°  =  e't 

2 


(2.4) 


A  restoring  couple  is  provided  by  the  tension  in  the  two  active 
sections ,  which  make  angles  0+0'  with  respect  to  the  center 
section.  For  small  0  and  6’,  the  restoring  torque  is 


t 


-2T 


sin  (Of  0*) 


=  -Tto  (ot  e') 


/  v 


(2.5) 


In  the  last  expression,  equation  (2.4)  has  been  used  to  elimi¬ 
nate  0'.  Thus  the  equation  of  motion  for  a  twisting  mode  is 


•  • 

10  *  T 

r 


(2.6) 


or 


u  20 
r 


=  o  »  where 


u 

r 


IAx 

i 


Here  1^  denotes  the  moment  of  inertia  appropriate  for  the  parti¬ 
cular  mode.  The  frequency  of  oscillation  is 


(2.7) 
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Like  the  frequency  of  the  translational  modes,  it  increases  as  the 
square  root  of  T. 

Typically,  the  parameters  of  the  pendulum- fiber  system  have 
the  following  values: 

m  =  35  8  t  =  2.5  cm 

T  =  lxlO6  dynes  £0=  7  cm 

From  equation  (2.3)  the  frequency  of  the  translational  modes 
is  estimated  as  24  Hz.  The  frequency  of  each  twisting  mode  depends 
on  the  moment  of  inertia  involved.  About  the  vertical  axis  of 
symmetry,  the  moment  is  1^  ««35g-cm2,  so  the  frequency  is  roughly 
110  Hz,  according  to  equation  (2  7).  The  moment  of  inertia 
about  the  horizontal  axis  perpendicular  to  the  fiber  at  its  center 
*2  =  I'^xlO^g-cm^.so  that  the  estimated  frequency  of  twisting 
about  this  axis  is  5.6  Hz. 

Owing  to  the  relatively  large  mass  and  length  of  the 
pendulum,  the  second  twisting  mode  has  quite  a  low  frequenev. 
Observations  of  the  fiber  pendulum-system  in  air  indicate  a  fre¬ 
quency  lower  than  the  calculated  value,  about  3  Hz.  Crude 
estimates  for  the  frequencies  of  the  translational  modes  from 
ccuation  (2.3)  are  also  in  error.  The  frequency  of  the 
horizontal  translation  has  been  measured  as  42  Hz  in  the  completed 
instrument.  The  first  twisting  mode  has  never  been  observed. 
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Although  in  the  vertical  accelerometer  one  of  the  unwanted 
modes  couples  directly  to  the  fundamental,6  this  is  not  the 
case  in  the  horizontal  instrument.  In  first  order,  the  unwanted 
modes  are  relatively  innocuous  in  the  horizontal  instrument. 

However,  non-linear  processes  i  the  fiber  can  still  cause  coupling 
of  energy  when  the  unwanted  modes  are  highly  excited.  This  is 
particularly  true  for  the  second  twisting  mode,  due  to  its  low 
frequency,  for  this  reason,  magnetic  damping  is  used  to  control 
the  excitation  of  the  unwanted  modes. 

The  apparatus  for  damping  the  troublesome  second  twisting 
mode  is  shown  in  Tigure  3.  Pendulum  damping  magnets  B13  and 

B14  produce  a  magnetic  field  which  passes  perpendicularly  through 
copper  pendulum  mass  B12.  Motion  perpendicular  to  the  field  lines 
causes  eddy  currents  which  dissipate  energy  by  ohmic  heating  in 
the  mass,  while  motion  along  the  field  lines  is  unaffected  in  first 
order.  Hence  the  twisting  mode,  but  not  the  fundamental  mode,  is 
damped.  The  vertical  translation  mode  is  also  damped.  In  order  to 
maximize  eddy  currents,  the  mass  is  made  of  copper.  The  magnets 
used,  which  have  six  pound  pull,  are  mounted  in  a  cylindrical  sleeve 
not  shown  in  the  figure. 

Preferential  damping  of  the  twisting  mode  depends  critical¬ 
ly  on  the  shape  and  thickness  of  the  mass  and  the  positions  of  the 
magnets.  Although  the  field  of  the  magnets  is  primarily  perpendicu¬ 
lar  to  the  pole  faces,  fringe  fields  are  appreciable.  Use  of  a 
thin,  flat  mass  in  the  region  of  the  most  linear  field  minimizes 
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the  effect  of  the  fringe  fields  in  damping  the  fundamental  mode. 
Spacing  of  the  magnets  is  determined  by  a  compromise  between  the 
strength  of  the  linear  field  at  the  center  and  the  volume  of  copper 
affected  by  fringing  fields.  Optimum  damping  is  obtained  when  a 
flat  mass  o.l  inch  thick  is  used  and  the  separation  between  the 
pole  faces  and  the  surface  of  the  mass  is  0.05  in. 

Other  magnetic  damping  assemblies  similar  to  those  used  in 
the  vertical  accelerometer  aid  in  the  elimination  of  unwanted  modes. 
Two  magnets  similar  to  B13  and  B14  are  mounted  on  either  side  of 
the  position  sensor  so  that  their  fields  pass  perpendicularly 
through  detector  plate  Bl.  These  magnets,  which  are  shown  in 
Figure  2  ,  t _  have  the  same  effect  as  the  pendulum  damping  magnets 
they  damp  the  vertical  translation  and  second  twisting  modes. 

With  the  two  systems  working  together,  the  Q  of  the  twisting  mode 
is  reduced  to  less  than  5,  while  the  Q  of  the  fundamental  mode 
is  20.  Auxiliary  damping  paddles  B8  and  BIO,  together  with  sixteen 
sr;all  magnets  not  shown,  reduce  excitation  of  the  horizontal  trans¬ 
lation  mode  and  the  first  twisting  mode.  With  this  auxiliary 
damping  assembly,  the  Q  of  the  horizontal  translation  is  about  30. 


3.  The  Position  Senior 


The  position  sensor  used  to  detect  motion  of  the  pen¬ 
dulum  is  identical  to  that  in  the  vertical  accelerometer.  It 
is  described  in  detail  elsewhere.7  The  essential  component  is 
a  three-element  parallel  plate  capacitor  whose  center  element 
is  aluminum  detector  plate  B1  in  Figure  3.  The  stainless  steel 
outer  plates  of  this  capacitor  are  stationary,  being  separated 
by  a  stable  quartz  spacer.  They  are  shown  in  Figure  2. 

The  outer  plates  of  the  sensing  capacitor  are  driven 
out  of  phase  at  a  frequency  of  16kHz  by  precisely  matched 
secondaries  of  a  toroidal  transformer.  The  resulting  signal 
on  the  detector  plate  is  proportional  in  amplitude  to  the 
displacement  of  the  center  plate  from  the  null  position.  After 
preamplification,  the  center  plate  signal  is  synchronously 
detected  with  respect  to  the  original  drive  signal.  The  output 
of  the  phase  sensitive  detector  is  then  a  measure  of  the  posi¬ 
tion  of  the  detector  plate.  As  a  whole,  the  position  sensor 
is  quite  linear  if  fringing  fields  in  the  capacitor  are  neglec¬ 
ted.  12  Its  sensitivity  is  high  enough  to  permit  detection  of  .05A  displace¬ 
ments  over  averaging  times  of  the  order  of  1  sec.7 

Placement  of  the  detector  plate  is  dictated  by  the 
requirements  of  ease  of  design.  Although  the  gain  of  the  sensing 
system  is  obviously  greater  if  the  sensor  center  plate  is  the 
pendulum  mass  itself,  the  sensitivity  of  the  sensor  is  sufficient 
to  allow  detection  at  the  present  position  without  serious  loos  of 
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gain.  The  advantage  accruing  from  such  a  design  is  that  the 
proven  structure  of  the  sensor,  fiber,  clamps  and  frame  from 
the  vertical  accelerometer  are  employed  without  alteration.  The 
gain  loss  is 

recovered  simply  by  increasing  the  electronic  amplification.  The  end 
clamps  and  tensioning  mechanism  for  the  fiber,  the  position  sensor  plates 
and  their  mounting  and  adjusting  apparatus,  and  the  supporting  frames  are 
the  same  as  in  the  vertical  instrument.  There  is,  however,  no  need  for 
external  adjustment  of  the  sensor  null  position,  as  the  horizontal  instru¬ 
ment  can  be  zeroed  by  tilting.  Thus,  the  parallel  motion  springs  upon 
which  the  outer  sensor  plate  assembly  ride;6are  locked  under  tension  at  the 
time  of  construction  and  no  mechanical  coupling  is  necessary.  The  assembled 
accelerometer  is  shown  in  Figure  5. 

Though  it  is  possible  to  apply  electrostatic  feedback  to  the 
detectoi  plate  to  close  the  accelerometer  loop,  all  the  preliminary 
data  have  been  taken  in  the  open  loop  configuration.  Feedback  opera¬ 
tion  is  accomplished  in  exactly  the  same  way  as  in  the  vertical  in¬ 
strument.7  Inverted  and  non-inverted  outputs  from  the  lock-in- 
amplifier  arc  applied  to  the  outer  plates  of  the  position  sensor.  In 
conjunction  with  a  small  bias  voltage  on  the  detector  plate,  these 
equal  and  opposite  feedback  voltages  produce  a  linear  electrostatic 
force  on  the  center  plate.  A  feedback  network,  consisting  of  an  in¬ 
tegrator  and  compensator,  tailors  the  gain  and  phase  of  the  loop  to 
make  the  system  stable.  The  integrator  provides  high  loop  gain  for 
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FIGURE  5 


Photograph  of  assembled  horizontal  accelerometer,  showing  mating 
frames  for  fiber-pendulum  assembly  and  outer  transducer  plates. 

Fiber  tensioning  plate  of  Figure  1  is  facing  camera. 

Insulated  plink  stoppers  in  center  of  outer  transducer  plates  are 
visible  inside  mode  damping  magnets.  Quartz  spacer  supporting 
outer  plates  is  attached  to  yoke  which  rides  on  parallel-motion 
springs  (not  visible)  at  right  and  left.  Since  transducer  zero  is 
adjustable  by  tilting,  bolt  in  bracket  under  thumb  at  left  locks  yoke 
in  position  under  tension  of  springs.  Supports  at  top  of  upper 
frame  attach  to  top  of  hard  vacuum  can.  At  right  is 
C.  William  Van  Sice,  who  assembles  instruments,  and  at  left  is  Jay 
Dratler,  Jr.  Barry  Block  is  taking  picture. 


stability  at  low  frequencies,  while  the  compensator  reduces  the 

resonance  peak  of  the  pendulum.  In  fedback  operation,  the  instrument 
output  is  taken  from  the  integrator. 

The  advantages  of  fedback  operation  are  stabilita. ion  of 
response  at  low  frequencies  and  reduction  of  actual  pendulum  motion, 
-hich  results  in  greater  linearity.  Because  the  fiber  is  torsion- 
ally  unstressed,  the  open-loop  linearity  of  the  horisontal  accelero¬ 
meter  is  a  priori  far  greater  than  that  of  the  open-loop  vertical 
instrument.  For  investigations  of  teleseisms,  long  term  stability 
Of  response  is  not  important,  so  feedback  is  of  little  value. 
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4.  Environmental  Control 


Like  its  vertical  counterpart, 6  the  horizontal  instru¬ 
ment  features  precise  internal  control  of  the  mechanical  system's 
environment.  Pressure  is  maintained  by  a  system  of  sealed  contain¬ 
ers,  while  temperature  is  regulated  by  a  fedback  temperature  control¬ 
ler  in  conjunction  witn  elaborate  thermal  insulation.  The  pendulum 
assembly,  fiber,  transducer,  and  supporting  frames  are  enclosed  in 

a  stainless  steel  hard  vacuum  can  sealed  with  a  gold  0-ring.  Pres- 

-7 

sure  in  the  can  is  10  torr.  During  evacuation,  the  can  and  its 
contents  are  baked  at  400°C  for  48  hours,  so  that  complete  :  ress 
relief  of  metal  parts  and  a  stable  high  vacuum  are  assured. 

Several  purposes  are  served  by  sealing  the  basic  accele,  - 
meter  in  a  high  vacuum  environment.  First,  isola¬ 
tion  from  external  pressure  fluctuations  is  obtained.  In  the 
absence  of  such  isolation,  changes  in  buoyant  forces  on  the  pendulum 
and  in  the  dielectric  of  the  position  sensor  constitute  sources  of 
noise.  Second,  the  mean  free  path  for  molecules  in  the  dilute  gas  is  made 
larger  than  a  typical  dimension  of  the  space  inside  the  contain¬ 
er.  Thus,  convection  and  its  attendant  problems  are  precluded. 

Third,  the  low  thermal  conductity  of  the  high  vacuum  provides 
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thermal  isolation  which  aids  in  stabilizing  the  temperature  of  the 
fiber. 

With  the  exception  of  part  of  the  temperature  control 
system,  the  hard  vacuum  can  and  its  contents  constitute  the  deploy¬ 
able  part  of  the  accelerometer.  The  small  diameter  (8  inches)  and 
light  weight  (20  lbs)  of  this  unit  make  it  suitable  for  deployment 
in  bore  holes  to  escape  the  effects  of  local  surface  tilting.  Such 
tilting,  which  is  the  most  formidable  source  of  earth  noise,  has 
been  observed  to  decrease  with  depth  of  deployment.1 9  For  bore 
hole  use,  redesign  of  the  temperature  control  system  for  compact¬ 
ness  is  necessary.  No  effort  has  yet  been  made  to  minimize  either 
the  size  or  weight  of  the  accelerometer  or  the  temperature  control¬ 
ler. 

For  deployment  in  vaults,  where  space  is  plentiful,  the 
environmental  control  is  completed  as  follows.  The  hard  vacuum 
can  is  enclosed  in  a  thick-walled  aluminum  soft  vacuum  can  sealed 
with  a  standard  rubber  0-ring.  Evacuated  to  a  pressure  of  1  torr, 
this  can  provides  additional  pressure  and  temperature  isolation. 

The  temperature  of  the  can  is  maintained  precisely  at  40°C  by  ohmic 
heating  in  a  bifilar  wire  wound  on  its  surface,  the  current  in 
which  is  supplied  by  a  remote  power  amplifier.  A  bifilar  winding 
is  used  to  eliminate  inductive  effects.  Heater  current  is  con¬ 
trolled  by  the  error  signal  from  a  precision  thermistor  bridge  and 
high  gain  preamplifier.  The  thermistor  is  mounted  directly  on  the 
heater  winding  with  thermally  conducting  epoxy.  Phase  compensation 
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and  response-tailoring  feedback  at  the  power  amplifier  insure  that 
the  thermistor  bridge,  preamplifier,  power  amplifier  and  heater 
form  a  stable  closed  loop.7  For  further  thermal  isolation,  the  soft 
vacuum  can  with  its  heater  winding  is  surrounded  by  several  inches 
of  polystyrene  balls  in  a  thermal  insulation  can.  Long  lightweight 
cables  (RG/187/U  coaxial  cable  and  30AWG  ribbon  cable),  coiled  in 
the  thermal  insulation,  provide  connections  between  the  soft  vacuum 
can  and  the  cables  leading  out  of  the  vault  while  maintaining  a  high 
resistence  to  heat  flow.  The  soft  vacuum  can  and  thermal 
insulation  can  are  mounted  on  a  remotely  controlled  tilting 
Platform  used  to  level  and  calibrate  the  instrument.  Long  legs  made 
from  thin-walled  stainless  steel  tubing  support  the  soft  vacuum  can 
without  ruining  the  thermal  isolation.  The  soft  vacuum  can,  thermal 
insulation,  temperature  controller,  and  tilting  platfonn  are  the 
same  as  in  the  vertical  instrument  The  only  difference  is  that  the 
soft  vacuum  can  in  the  horizontal  has  longer  legs  and  a  snoot  to 
house  the  protuberance  on  the  hard  vacuum  can. 

5-  Temperature  Stabilization 

In  the  vertical  accelerometer,  precise  temperature  control 
is  needed  because  the  torsion  in  the  fiber,  which  has  a  large  tempera¬ 
ture  coefficient,  balances  the  entire  force  of  gravity.  This  diffi¬ 
culty  does  not  occur  in  the  horizontal  instrument,  where  quiescent  tor¬ 
sion  in  the  fiber  is  small  (See  section  1).  Theoretically,  the 
horizontal  accelerometer  has  zero  temperature  coefficient.  In 
practice,  however,  a  temperature  dependence  arises  from  errors  in 
setting  the  fiber  at  its  torsionally  unstressed  position  during 
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assembly.  Due  to  such  errors,  there  is  some  small  torsion  in  the 
fiber  when  the  instrument  is  level.  Changes  in  this  torsion  with 
temperature  produce  spurious  output. 

A  worst-case  analysis  proceeds  as  follows.  The  maximum 
error  in  installing  the  fiber  is  0g  *  0.1  radian  of  rotation. 

The  tortion  constant  for  the  fiber  in  the  vertical  accelerometer 
is  given  in  Section  1.  Since  the  initial  rotation  of  the  fiber  in 
the  vertical  is  considerable,  the  torsion  constant  for  small  addi¬ 
tional  rotations  must  be  larger  than  the  torsion  constant  of  the 
nearly  unstressed  fiber  in  the  horizontal.  Thus,  a  reasonable 
upper  bound  for  the  torsion  constant  in  the  horizontal  is 
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7. 5x10 3  dyne -cm 

rad 


(4.1) 


The  corresponding  upper  limit  for  the  quiescent  fiber  torsion  in 
the  levelled  horizontal  instrument  is  =  7.5xl02  dyne-cm. 

Using  the  values  of  the  mass  m  and  effective  torsional  lever  arm 
Lm  of  the  pendulum  assembly  given  in  equations  (1.4),  the  equiva¬ 

lent  horizontal' acceleration  is  estimated  as 
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1.4cm/sec2 


10-3g 


(4.2) 


The  fiber  positioning  error  in  the  horizontal  produces  a  maximum 
.  .  -3 

effective  acceleration  of  10  g.  The  temperature  coefficient 
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control  ,,  o^.*  for  occurntr  trophy,  I  col  „„orw„„. 

T.  Mun  t.lMtlMlc  UMlmtloM  of  10*“|,for  Inst.nec,  the 
fiber  tMp#retur#  euet  he  stable  a,  th*  niil {degree  level.  This 
temperature  regulation  is  aeconplished  in  tvo  ways.  rirst,  the 
eechamcal  systee  Is  isolated  fron  environmental  temperature  changes 
hy  th#  systen  of  passive  thermal  Insulation  consisting  of  the  hard 
and  soft  vacuum  cans  and  the  polystyrene  bails.  Second,  the  temper¬ 
ature  of  the  control  surface  surrounding  the  soft  vacuun  can  is  held 
fixed  by  the  fedback  temperature  controller.  The  vacuum  cans  provide 
th«mal  filtering  between  the  fiber  and  the  control  surface,  while 
the  polystyrene  balls  isolate  the  control  surface  from  the  outside  world. 


In  coftftUfrlnf  tU  i  to  mail  IfwjiaUcn,  An  analogy  with 
f  l*jwp**A  filters  in  alactronics  it  useful.  ?x  §yfct*m  with  a 
|jlr**  Mai  cecity  is  contAln^  in  *  can  with  lew  itomal  e  m- 
ihict iv*ty,  i.».,  high  ttomal  resistance.  T**f*er«uure  fluciwa- 
tion*  in  the  external  environment  create  a  teeperature  differ¬ 
ence  between  the  Inside  end  outside  of  the  cen.  This  difference 
drive®  a  t hemal  current,  rwneiy  a  heat  flow,  into  the  inner  sys- 
te«.  The  heat  flow  is  related  to  the  temperature  difference 
by  the  integrated  thereat  resistance  of  the  outer  can,  which 

is  analogous  to  the  resistance  in  Ohm's  Uw.  As  the  thermal  current 
heats  the  inner  system,  it  causes  a  te-p« rat or* 

rise  proport i >nal  to  its  integral  over  tire  and  inversely  propor¬ 
tional  to  the  system's  heat  capacity.  This  is  analogous  to  the 
charging  of  a  capacitor.  Taken  as  a  whole,  the  thermal 

system  is  like  an  RC  loupass  filter.  When  a  temperature  step  is 
generated  outside  the  can,  some  time  elapses  before  enough  heat  can 
flow  through  the  can's  resistance  to  "charge"  ,he  lnMr  syEte„  to 
a  higher  temperature.  The  time  required  for  1/e  of  the  external 
temperature  step  to  appear  at  the  inner  system  is  the  thermal  time 
constante  It  increases  directly  as  the  thermal  resistance  of  the 
container  and  the  thermal  capacitance  of  the  inner  system.  As  in 
the  electronic  lowpass  filter,  any  sinusoidal  fluctuation  in 
external  temperature  whose  period  is  appreciably 

shorter  than  the  time  constant  is  considerably  attenuated  at  the 
inner  system. 
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4!M)' ty  Utv««n  ichvnMl  and  electronic  filters,  however, 
4n,m  b#  art  ed  too  for.  Tho  p4ra««(eni  of  ih*  electronic  circuit 
«r*  localised  In  dlfcret*  component*,  while  the  thermal  resistance 
«r  4  opoci  one#  or*  mixed  and  distributed  over  space.  The  outer 
ontalner,  which  acts  an  tha  thermal  resistor,  has  a  fair  amour,: 
of  heat  capacity,  and  tha  thermal  capacitor  has  s^me  thermal 
resistance.  Both  the  thermal  resistance  and  capacitance  depend  on 
th*  detaIls  of  construction,  materials,  and  parts  placement .  Tor 
this  reason,  such  thermal  systems  are  quite  difficult  to  analyze 
mathematically.  The  general  idea  of  a  lovpa3c  filter  with  a  charac¬ 
teristic  time  constant  is  useful  in  estimating  the  effectiveness  of 
thermal  insulation. 

In  the  quartz  fiber  accelerometer,  there  are  three  stages 
of  thermal  filtering.  The  high  ti.jrmal  resistance  of  the  polystyrene 
balls,  followed  by  the  large  heat  capacity  of  the  soft  vacuum  can, 
constitutes  one  stage  of  filtering.  The  next  stage  relies  on  the 
thermal  resistance  of  the  soft  vacuum  and  the  thermal  capacity  of 
the  hard  vaccum  can.  Finally,  the  thermal  resistance  of  the  hard 
vacuum,  coupled  with  the  thermal  capacity  of  the  mechanical  system 
and  its  frame,  provide  the  third  stage  of  filtering. 

Reductions  in  the  effective  thermal  resistances  due  to 
alternate  paths  of  heat  flow  -'e  quite  important.  The  polystyrene 
insulation,  for  instance,  is  bypassed  by  a  thermal  path  through 
the  tilting  platform  and  up  the  three  legs  which  support  the  soft 
vacuum  can.  These  legs  are  made  of  thin-walled  stainless  steel 
tubing  to  maximize  the  thermal  resistance  of  this  path.  Similarly, 
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it  vacuum  can  Inert...  ,h.  thermal  rt.i.tanc.  of  this  path.  The 

hard  vacuum  1.  bypaseed  by  th.  short  la*,  with  which  tho  mechanical 
system  and  it,  frame  are  suspended.  Since  these  legs  are  not 
insulated,  the  final  stage  of  thermal  filtering  is  not  as 
effective  as  the  others.  Another  important  thermal  path  is 
through  the  tilting  platform,  which  is  a  large,  thick  plate  of 
aluminum.  To  protect  it  from  direct  exposure  to  the  thermal  environ 
ment,  it  is  usually  surrounded  by  polystyrene  balls  contained  in 
a  circular  layer  of  foam  rubber  on  the  vault  floor. 

Perhaps  the  most  damaging  alternate  path  of  heat  flow  is 
through  the  electrical  cables  to  the  mechanical  system  and  tempera¬ 
ture  sensor.  In  order  to  maximi.e  the  effective  thermal  resistance 

these  cables,  lengths  of  them  are  coiled  inside  the  thermal 
insulation  can. 

In 

addition,  lightweight  cables  are  used  to  connect  the  mechanical 
system  to  the  external  electronics.  The  coaxial  cables  are  teflon- 
insulated  RG187/U  ,  while  the  other  conductors  are  part  of  a 

30AWG  ribbon  cable.  After  emerging  from  the  thermal  insulation 
can,  these  cables  are  connected  to  heavier  standard  cables  (RG62/U 
and  18A-.G  multiconductor),  which  connect  the  deployable  unit  to 
the  electronics  panel.  With  the  thermal  insulation 


constructed 

between  the 


as  outlined  ,  the  time 
outside  world  and  the  fiber 


constant  for  heat  transfer 
is  several  hours.  Thus, 
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tho  thermal  filtering  drastically  reduces  the  effects  of  external 
temperature  fluctuations  with  periods  of  a  few  minutes  or  less. 

Lower  frequency  variations  are  eliminated  by  the  temperature  controller. 


Figure  6  shows  a  flow  chart  of  the  fedback  tempera¬ 
ture  controller.  The 

thermistor  measures  the  temperature  at  a  point  on  the  control  sur¬ 
face.  The  temperature  sensor  and  preamplifier  compares  the 
termistor  resistance  to  that  of  a  standard.  When  the  two  differ, 
it  produces  an  error  signal  proportional  to  the  difference.  This 
temperature  error  signal  is  sent  to  the  fedback  power  amplifier, 
whose  outpu.  determines  the  current  in  the  resistive  heater 
wound  around  the  control  surface.  The  error  signal  changes  the 
heater  current  to  bring  the  control  surface  temperature  back  to 
the  null  point.  Since  the  null  point,  which  is  determined  by 
the  resistance  of  the  standard  in  the  sensor,  exceeds  room 
temperature,  there  is  always  current  in  the  heater  winding.  This 
mean  current  compensates  for  thermal  leakage  through  the  insula¬ 
tion.  Fluctuations  in  the  heater  current  maintain  the  thermistor 

temperature  fixed  within  4  x  10  ^°C  per  100  sec 

*"8 

and  4  x  10  °C  per  day.7 
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FIGURE  6 


Flow  chart  of  fedback  temperature  controller  for  accelerometer. 
Control  surface  is  metallic  cylinder  surrounding  mechanical 
instrument.  It  is  thermally  insulated  from  both  fiber  and 
external  environment. 
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In  practice,  tha  active  temperature  controller  and  the  passive 
thermal  insulation  interact  considerably.  Not  only  does  the  in¬ 
sulation  determine  the  mean  heater  current  needed  to  maintain  a 
given  temperature,  but  it  also  affects  the  stability  of  the  tempera¬ 
ture  controller .  For  small  fluctuations  in  the  heater  current ,  it 
controls  the  loss  of  heat  from  the  surf-.ee  and  thus 

limits  the  thermal  lag  between  the  heater  and  thermistor. 

Since  this  thermal  lag  is  primarily  responsible 

for  difficulties  in  stabilizing  the  fedback  controller,  improvement 
in  the  thermal  insulation  is  thus  one  of  the  expedients  to  be  used 
in  stabilizing  recalcitrant  control  systems.  As  in  the  case  of  the 
design  of  the  insulation  itself,  calculations  of  the  stabilizing  effects 
are  made  prohibitively  difficult  by  the  fact  that  the  parameters 
are  distributed.  Such  improvements  an  most  quickly  made  by  trial 
and  error. 

The  temperature  controller  can  increase  the  need  for  thermal 
insulation  by  providing  small  thermal  fluctuations  of  its  own. 

As  heating  depends  on  the  square  of  the  current,  the  controller  as  a 
whole  is  a  non-linear  system.  It  is  subject  to  marginal  instabilities. 

Usually,  these  are  at  rather  high  frequencies  where  the  controller 
loop  gain  is  low.  In  the  presence  of  strong  environmental  tempera¬ 
ture  drift,  however,  the  controller  sometimes  cycles  at  the 
higher  teleseisrr.ic  frequencies.  To  reduce  the  effect  of  these 
low-level  thermal  oscillations,  the  control  surface  is  placed  on 
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tho  outside  of  the  soft  vacuus  can,  where  it  is  better  insulated  free 
the  fiber. 

When  tho  temperature  controller  is  working  properly,  it  com¬ 
plements  the  thermal  insulation  by  providing  regulation  for  periods 
longer  than  the  time  constant  of  the  .nermal  filtering.  The  loop 
gain  of  the  temperature  controller  is  quite  high  below  10’4  Hz  and 
falls  20  db/dec  above  this  frequency. 

The  active  device  thus  provides  good  regulation  where  the  passive 
insulation  is  ineffective.  Conversely,  the  insulation  attenuates 
temperature  fluctuations  at  higher  frequencies,  where  the  low  loop 
gain  reduces  the  effectiveness  of  the  temperature  controller.  The 
frequency  region  of  worst  temperature  regulation  is  in  the  neighbor¬ 
hood  of  100  sec  periods,  where  neither  the  active  nor  the  passive 
system  is  optimally  effective. 

The  overall  effectiveness  of  the  temperature  regulation 

is  best  measured  by  the  quality  of  the  accelerometer  output.  In 

the  vertical  accelerometer 
the  seismic  spectral  region,  data  from  /  indicate  that  the  total 

instrumental  noise  level  is  below  ground  noise  levels,  which  are 

of  the  order  of  10  g.  Thus  stabilization  of  the  fiber  temperature 

over  teleseismic  periods  (20-50  sec)  is  good  to  a  microdegree, 
vertical 

As  for /tidal  observations,  they  agree  with  theory  to  better 
than  1%.  This  implies  fiber  temperature  control  of  the  order 
of  10  °C  over  half-day  periods.  The  only  exception  to  the  generally 
high  degree  of  stability  is  in  the  neighborhood  of  100  sec  periods. 
Some  unexplained  noise  has  been  observed  in  this  region.  However, 
the  origin  of  this  noise  is  not  known,  since  ground  noise  levels 
at  these  frequencies  have  not  been  accurately  determined. 
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L  ,  the  torque  is 

Ti  e  tMhlmcos0  s  ^m\  (6.n 

The  restoring  torque  due  to  gravity  is  given  in  equation  (1.1),  The 
damping  torque  is  simply  given  by 

Td  =  "b®  (6.2) 
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where  To  •  —  in  the  fr*e  period. 

This  gain  has  the  same  form  as  t*v*t  for  the  v«rtiaal 

mass-fiber  system.  Inc*  L  is  *4de  th*  troe  In  ioth  t«firia<«f.3i. 

the  horisontal-to-vert leal  OC  gain  rails  U  given  tty 


(6.7) 


where  the  subscripts  v  and  h  denote  the  obvious. 

Numerical  values  of  the  quantities  in  espresticn  (4.7) 
are  obtained  by  a  combination  of  calculation  and  eeasureeent .  The 
masses  are  determined  by  direct  weiring  of  parts,  while  the  »o- 
ments  of  inertia  are  calculated  from  the  known  weights 

and  dimensions.  Moment  arm  is  determined  by  balancing  the 
pendulum  assembly,  and  Lav  is  calculated  free  part  weights  and 
blueprint  dimensions.  The  free  periods  are  e^a  ured  ,froa  chart 
recordings  of  the  decaying  fundamental  sode  oscillations.  In  Table  1, 
parameter  values  for  a  typical  vertical  instrument  and  the  prototype 
horizontal  instrument  are  listed.  The  gain  ratio  (5.7  )  is  computed  as 
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The  DC  open  loop  gains  have  been  measured  for  both  instruments.  Their 
ratio  is  0.13,  which  is  26%  below  the  value  calculated  from  crude 
measurements  o.  the  moments 

Th*re  is  one  important  difference  between  the  vertical  and 
horizontal  instruments  which  is  not  illuminated  by  response  calcu¬ 
lations.  Below  the  resonance  frequency,  the  horizontal  is  a  tilt- 
meter,  while  the  vertical  measures  tilt  only  in  second  order.  In 
deriving  the  response  6.5  for  the  horizontal,  it  was  assumed 
that  the  instrumental  and  local  gravitational  verticals  were  coin¬ 
cident.  If  this  is  not  true,  there  may  be  a  non-zero  component 
of  the  acceleration  of  gravity  along  the  degree  of  freedom  of  the 
accelerometer.  Such  a  component  of  gravity  along  the  instrumental 
horizontal  depends  on  the  sine  o'  the  tilt  angle,  so  there  is  a 
linear  relationship  between  tilt  and  accelerometer  output  for 
small  DC  tilts.  The  output  of  the  horizontal  accelerometer  can  thus 
be  calibrated  in  terms  of  volts/g  of  acceleration  or  volts/radian  of 
tilt.  For  small  tilts,  these  units  are  equivalent,  as  the  effects 
on  the  instrument  of  tilt  and  kinematic  horizontal  acceleration  are 
physically  indistinguishable.  As  for  the  vertical  instrument,  the 
component  of  gravitational  acceleration  along  its  degree  of  freedom 
depends  on  the  cosine  of  the  tilt  angle,  so  that  it  responds  in 
second  order  to  small  tilts.  The  DC  responses  of  the  two  instruments 
to  tilt 

are  used  in  calibration.7  Above  their  resonance  frequencies,  both 
instruments  respond  in  first  order  to  tilt  because  the  passive 
elements  (mass  or  pendulum)  cannot  respond  with  sufficient  speed, 
so  they  remain  virtually  stationary  while  the  outer  sensor  plates 
oscillate  about  them. 
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7*  Filtering  and  Recording  of  Data 

For  geophysical  purposes,  the  quartz  fiber  accelerometers 
have  broadband  response.  In  the  open  loop  system,  the  region  of 
high  gam  extends  from  DC  to  the  resonance  frequency. 

The  open  loop  bandwidth  is  thus  about  1  Hz.  Such  band¬ 
width  is  sufficient  for  simultaneous  investigations  of  earth  tides, 
free  oscillations  of  the  earth,  and  long  period  teleseisms.  Very 
often,  however,  it  is  not  convenient  to  use  the  same 
data  record  to  study  all  three  phenomena.  Tidal  data  appears  cleaner, 
and  is  more  easily  analyzed,  when  earthquakes  are  removed  from  it. 
Study  of  teleseismic  records  is  only  hampered  by  apparent  drift  due 
to  earth  tides.  Certain  studies  require  higher  gain 

than  that  available  at  the  accelerometer  output.  For  these  reasons, 
the  raw  output  of  the  accelerometer  is  modified  by  three  active 
analog  filters,  which  respond  preferentially  to  tides,  free  oscilla¬ 
tions,  and  teleseisms.  Because  they  modify  the  accelerometer  output 
before  it  is  recorded,  these  devices  are  called  prefilters. 

Each  of  the  three  prefilters  is  an  electronic  circuit  which  uses 
one  or  more  fedback  operational  amplifiers  to  set  its  response.7  Each 
has  a  full  scale  output  of  ±10V.  The  tidal  prefilter  simply  attenuates 
high  frequency  signals,  while  the  normal  mode  and  seismic  prefilters 
both  have  gainful  response  at  certain  frequencies. 
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The  tidal  prefilter  is  a  double-pole  lowpass  filter.  Its  cut¬ 
off  frequency  can  be  made  5.75  cph  (slow)  or  57.5  cph  (fast)  by  means 
of  a  switch,  for  tidal  investigations,  the  lower  cutoff  frequency 
is  used.  Since  the  highest  tidal  frequency  of  usual  interest  is 
that  of  the  triurnal  tide,  or  0.175  cph,  the  slow  cutoff  frequency 

gives  the  tidal  prefilter  nearly  unity  gain  and  tore  phase  shift 
over  the  entire  tidal  band. 

In  calibration  and  testing  of  the  instrument,  the  response 
near  DC  is  of  major  interest.  The  tidal  prefilter  is  used  during 
calibration  to  observe  DC  levels  while  averaging  out  fluctuations 
Cue  to  microseismic  ground  „oise.>°To  avoid  long  waiting  periods, 
the  time  constant  is  set  to  fast  for  such  work. 

For  studying  the  normal  modes  of  free  oscillation  of  the 

earth,  a  gainful  bandpass  filter  is  uoH  t«- 

P  filter  is  used.  Its  passband  is  variable, 

having  ,  lower  knee  frequency  of  1  cph  and  an  upper  knee  frequency 

Of  30.60.  120  or  180  cph.  The  gain  th.  passband  i,  UOdb  above 

the  DC  level.  On  the  high  side  of  the  passband.  the  response  fall, 
UOdb/dec.  On  the  low  side,  it  falls  UOdb/dec  to  0.1  cph  and  is 
constant  thereafter.  g„  input  attenuator  sets  the  DC  level.  When 
the  attenuation  is  zero,  the  gain  at  DC  is  unity. 
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The  normal  mode  prefilter  is  usually  run 
with  no  input  attenuation  and  the  upper  knee  frequency  set  at 
30  cph.  This  configuration  provides  sucficient  gain  and  bandpass 
for  studies  of  the  lower  frequency  normal  modes.  Increasing  the 
upper  knee  frequency  includes  more  of  the  normal  mode  band  in  the 
gainful  region.  At  the  highest  setting,  long  period  seismic  signals 
are  included  in  the  bandpass.  This  setting  allows  all  seismic 
signals  with  periods  greater  than  20  sec  to  be  observed  with  nearly 
the  same  gain.  Unfortunately,  large  Rayleigh  waves  and  ground  noise 
tend  to  saturate  the  prefilter  electronics  when  its  bandwidth  is 
this  large.  This  is  the  reason  the  input  attenuator  is  provided. 

It  allows  a  compromise  between  gain  and  bandwidth  to  be  made.  The 
nature  of  this  compromise  depends  on  the  size  of  events  to  be 
observed  and  the  level  of  ground  noise  at  the  instrument  site. 

The  reason  for  the  relative  attenuation  below  0.1  cph  is  the 
large  amplitude  of  the  tidal  signals.  Peak-to-peak  tidal  accelera¬ 
tions  can  be  as  high  as  3xl0~7g,  while  normal  modes  and  seismic 
waves  typically  produce  accelerations  of  order  10  ^g.  Tidal  signals 
at  the  lock-in  amplifier  output  are  thus  near  IV  peak-to-peak.  If 
amplified  significantly,  these  signals  saturate  the  prefilter. 

Phase  shifts  due  to  the  changing  gain  of  the  pre filter  near  tidal 
frequencies  displace  and  distort  the  tidal  signals.  This  precludes 
use  of  the  normal  mode  prefilter  output  for  tidal  studies  or  tidal 


calibration. 


The  seismic  prefilter  is  a  peaked  bandpass  filter.  Its  peak 

is  fixed  at  100  second  periods,  while  its  half  power  points  are  at 

60  and  200  second  periods.  At  high  frequencies,  its  response  falls 
n 

as  hi  ,  or  140db/dec,  while  at  low  frequencies  its  response  dies 
2 

as  u  ,  or  40db/dec.  The  peak  gain  is  set  by  the  input  attenuator, 
which  has  six  positions,  and  can  be  varied  from  60  to  5830.  Normally, 
the  attenuator  is  set  for  a  peak  gain  of  621. 

The  high  peak  gain  and  the  rather  narrow  bandwidth  of  the 
seismic  prefilter  are  quite  useful  in  observing  long  period  tele¬ 
seisms.  They  permit  measurement  of  seismic  signals  as  low  as  10  ^°g 
in  amplitude.  The  sharp  rolloff  of  the  response  at  high  frequencies 
reduces  interference  due  to  tnicroscistnic  signals,10  while  the 
rolloff  at  low  frequencies  eliminates  the  tides.  The  seismic  pre¬ 
filter  is  used  primarily  for  detection  of  distant  teleseismic  events. 

The  accelerometers  have  flat  responses  to  acceleration 
over  the  bandpasses  of  the  prefilters.  Thus,  the  shape  of  the 
composite  response  at  the  output  of  each  prefilter  is  determined 
by  the  response  of  the  prefilter  itself.  The 

gain  level  is  determined  by  the  DC  gain  of  the  accelerometer.  The 
exact  form  of  the  response  to  acceleration  at  each  prefilter  output 
is  calculated  by  multiplying  the  prefilter  response  function  by  the 
appropriate  accelerometer  response.  The  resultant  total  response 
functions  are  plotted  along  with  the  accelerometer  response  for 
the  open  loop  system  in  Figure  7. 
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FIGURE  7 


Composite  responses  to  acceleration  of  open  loop  instrument 
and  three  prefilters,  with  response  at  position  sensor  output 
shown  for  comparison.  All  response  magnitudes  are  plotted  relative 
to  open  loop  DC  gain  of  accelerometer. 


FREQUENCY  (Hi! 


The  output  of  each  prefilter  is  recorded  digitally  on  magnetic 
tape  with  a  resolution  of  +16384  digitizer  counts  per  +10  volts. 

In  the  cases  of  the  normal  mode  and  seismic  prefilters,  the  outputs 
with  +10V  full  scale, 

are  used  for  digitizing.  Although  various  sampling  intervals 
are  available,  the  typical  interval  is  1  second,  so  that  the  digi¬ 
tizer  Nyquist  frequency  is  1/2  Hz.  Part  of  the  function  of  the 
prefilters  is  to  provide  plenty  of  attenuation  above  the  Nyquist 
frequency  to  eliminate  aliasing  of  energy  in  the  digitizing  process. 
As  Figure  7  shows,  the  response  of  each  prefilter  at  the 
Nyquist  frequency  is  down  at  least  80db  from  its  peak  response.  Such 
attenuation  is  sufficient  to  eliminate  aliasing  for  all  practical 
purposes . 

It  is  impossible  to  overemphasize  the  importance  of  digital 
recording  of  data.  Once  a  digital  tape  recording  is  made,  all  the 

computational  power  of  high  speed  electronic  computers  is  readily 
available.  This  puts  the  mathematical  tools  of  modern  data  analysis 
within  easy  reach.  In  particular,  digital  filtering  and  fast- 
Fourier  transformation  allow  quick  and  quantitative  determination 
of  the  spectral  content  of  various  signals.  Cross -correlation  aids 
in  establishing  a  relationship  of  cause  and  effect.  And  high- 
resolution  computer  plots  provide  a  simple  and  accurate  means  of 
displaying  results. 
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Unfortunately,  many  geophysicists  have  realited  too  late  that 
the  difference  between  digital  and  analog  records  is  not  merely 
quantitative.  The  dynamic  range  of  digital  records  is  as  high  as 
the  digitiser  resolution,  usually~lo\  Due  to  the  width  of  the 
pen  and  the  servo  dead  band,  the  dynamic  range  of  a  chart  recorder 
is  seldom  much  larger  than  200.  This  lower  resolution  allows  much 
of  the  information  in  small  signals  to  be  buried  in  the  pen  width. 
Once  lost  in  a  chart  record,  such  information  can  never  be  recovered. 
In  addition,  the  effort  and  expense  involved  in  digitising  even 
large  signals  from  analog  records  is  often  prohibitive. 

Chart  records  are  valuable  as  diagnostic  tools  and  catalogues 
of  events,  but  digital  records  are  indispensable  for  quantitative 
work.  Without  direct  digital  records,  the  analysis  presented  here 
and  elsewhere  could  not  have  been  attempted. 


.118 


52 


8.  Noise 


There  are  three  sources  of  noise  in  quartz-fiber  accelerometer 
systems:  Brownian  motion  of  the  mass  or  pendulum,  fluctuations  01 
the  fiber  temperature,  and  spurious  signals  in  the  sensor  and  pre- 
filter  electronics.  Reduction  in  noise  from  the  latter  two  sources 
is  largely  a  function  of  the  effort  expended  in  perfecting  the 
temperature  controller  and  electronic  systems.  For  all  accelerometers 
the  Browian  noise  level  represents  an  intrinsic  limit  to  signal 
detection.  It  depends  on  the  construction  and  operating  temperature 
of  the  mechanical  accelerometer,  but  is  otherwise  invariable. 

In  the  quartz-fiber  accelerometers.  Brownian  noise  is  pro¬ 
duced  by  collisions  of  gas  molecules  with  the  mass  or  pendulum, 
random  motions  of  electrons  in  the  damping  plates  through  the 
fringe  fields  of  the  damping  magnets,  and  fluctuations  in  the 
crystal  structure  of  the  mass  or  pendulum  assembly.  Because  these 
random  phenomena  are  uncorrelated  in  time,  their  effects  can  be 
expressed  as  an  equivalent  input  acceleration  with  a  flat  spectrum. 

The  constant  power  spectral  density  P  (w)  of  the  Brownian  acceleration  noise 
depends  on  the  mechanical  parameters  of  the  accelerometer  and  the 
absolute  temperature  of  the  immediate  environment.  The  formula  is^ 


P  (w) 
a 


—  ^  4KnT 

C  2  QI  B 
m 


(8.1) 
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where  l,  Ca,  «0,  Q,  «nd  I  are  *5  |n  equation  6.3  ana  &.w, 

Kp  1.38x1  1  ule/°K  it  toltssann's  constant,  and  i 

it  the  absolute  temperature* 

Except  for  the  focartrical  factor  l.f/C*  ,  formula  g.i 
applies  to  all  similar  accelerometers.  The  Brownian  noise  poutr 
varies  inversely  as  the  free  period  T0  -  ?»/w0,  the  mc^nt  of 
inertia  I,  and  the  quality  factor,  Q.  Aa  #hown  by  *qu4tlen  6 

the  DC  gain  of  the  acceleroaieter  depends  on  the  square  of  the 

free  period.  Thus,  a  large  free  period  provides  both  high  gain 
and  low  noise.  Tor  this  reason, 


geophysical  accelerometers  have  traditionally  been  de¬ 
signed  with  masses  that  can  only  be  described  as  gargantuan.*  In 
the  quartz  fiber  accelerometers,  the  free  period  Is  made  suffi¬ 
ciently  large  by  using  a  weak  restoring  device,  and  the  moment 
of  inertia  is  made  rather  small.  Thus  unwieldy  size  and  weight 
are  avoided.  The  Q  of  the  accelerometer,  which  is  still  left  as 
a  free  parameter,  is  made  large  enough  to  reduce  Brownian  noise 

to  an  acceptable  level.  In  practice,  the  Q  is  determined  by  fringe 
Helds  of  the  damping  magnets  and  not  by  the  resldual  ^ 

the  hard  vacuum  can.  If  operation  in  areas  with  low  site  noise  is 
contempleted ,  the  5  can  be  raised  by  increasing  the  dating  magnet 
spacing  and  thus  reducing  all  damping,  otherwise,  careful  work 
is  required  to  maximize  the  Q  of  the  fundamental  mod. 
while  minimizing  the  Q's  of  all  the  unwanted  modes.  Normally, 
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r  'h*  r4#r  of  ?0  «* I  th*  brownian  noif*  level  a*  about  lO’^g'/cph, 
* ritf*  U  **f **«*««!  ly  lew  for  **t  geophysical  observai loos. 

alculaiions  of  the  fcrewr.Un  note*  levels  nr*  readily 
****  trm  f  rmji4  ••*•  Th*  necessary  tnsiru**r,*al  parameters 
r°r  'iy?Ic4l  vertical  and  horlsontal  accelerometers  are  giver,  Jr, 
r*M#  i.  Uemally,  the  temperature  U  T  *  313*K.  The  power  spectral 

teMltles  of  brownian  acceleration  at  the  inputs  ar*  list*)  In 
Table  ?. 

TAE1X  ? 

Vertical  Moriaontal 

COJ)  (HD 

*a  *  g  /cph  p*  «  I.S*10’*^g^/cph 

•  T.MIO  ?lg?/Hs  •  3.0k10"7Ie^/Hz 

* r°*  th#e«  v*iu#»,  the  noise  levels  at  any  point  in  the 
accelerometer  electronics  due  to  Brownian  noise  «t  the  input  ma y 
b*  calculated.  Over  the  bandwidth  of  Hot  response,  the  power 
spectral  density  of  fluctuation*  at  the  position  sensor  output  are 
determined  by  multiplying  the  number*  in  Table  ?  by  the  square 
of  the  DC  accelerometer  gain.  This  gain,  which  is  taVen  as 
unity  in  rigure  7,  is  determined  by  static  lilt  calibration  of 
the  instrument.'’  Tor  typical  vertical  and  horizontal  accelero¬ 
meters  the  power  spectral  density  of  fluctuations  at  the  position 
sensor  output  arc  given  in  Table  3. 
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TABLE  3 


Vertical  Horizontal 

(VI. 3) 

(Q3)  (HI) 

Py(sens)=  8.2xl0"nv2/cph  Pv(sens)=  2.3xl0"12v2/cph 

=  2.9x10"7v2/Hz  =  8.0x10"9v2/Hz 

Of  course,  these  values  apply  only  below  the  resonance  peak  in  the 
open  loop  response,  i.t. t  in  the  region  of  flat  response. 

Since  the  p^filters  have  appreciable  gain  only  in  the 
region  of  flat  response,  the  accelerometer's  resonance  peak  and 
the  filtering  in  the  phase  sensitive  detector  can  be  ignored  for 
the  purpose  of  calculating  fluctuations  at  the  prefilter  outputs. 

Each  prefilter  is  considered  subject  to  white  noise  at  its  input. 

The  power  spectral  density  at  the  output  is  obtained  by  multi¬ 
plying  the  input  density  by  the  square  modulus  of  the  prefilter 
response  function.  The  output  spectrum  thus  assumes  the  shape 
of  the  squared  filter  response.  By  integrating  the  outPut  spectral 
density  over  frequency,  the  mean  square  value  of  voltage  fluctuations 
at  the  prefilter  output  is  obtained.  The  mean  square  and  RMS 
fluctuations  at  the  three  prefilter  outputs  are  given  in  Table  4 
for  both  the  vertical  and  horizontal  accelerometers. 
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Fluctuations  at  Prefilter  Outputs  Due  to  Brownian  Noise 
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In  addition  to  Brownian  noise  at  the  input,  there  are 
several  sources  of  spurious  signals  in  the  accelerometers.  The 
most  common  of  these  is  Nyquist  noise  in  resistors.  Also,  time- 

varying  magnetic  fields  in  the  environment  can  produce  large  and 
meaningless  signals  if  ground  loops  are  present.  These  arc  the 

standard  ills  to  which  all  electronic  equipment  is  prey.  Special 
evils  occur  in  electronics  designed  for  use  at  low  frequencies.  One 
of  these  is  1/f  noise,  which  is  common  in  DC  amplifiers.  At  the  extremely  low 
frequencies  involved  in  geophysical  observations,  temperature 
coefficients  of  components,  particularly  those  of  diode  and  tran¬ 
sistor  junctions,  are  another  source  of  noise. 

In  the  position  sensor  circuitry,  the  vicissitudes  of  low 

frequency  work  are  avoided  through  the  use  of  phase  sensitive 
detection.  The  preamplifier  and  leading  stages  of  the  lock-in¬ 
amplifier  operate  at  16  KHz  ,  so  that  only  Nyquist  noise  and 
ground  loops  are  significant.  Theoretically,  Nyquist  noise  in  the 
lOOMft  resistor  across  the  preamplifier  input 
constitutes  the  essential  limitation  of  the  sensor  electronics. 


According  to  the  Nyquist  formula,  the  power  spectral  density 


of  random  voltage  fluctuations  in  a  resistor  is  13. 


where  Kg  -  1.30x10  joule/°K  is  Boltzmann's  constant,  T  is  the 
absolute  temperature,  and  R  is  the  resistance.  At  the  controller 
temperature  of  313°K,  the  noise  power  density  produced  by  the  108 
ohm  resistance  is  about 


P 

v 


*  1.7xlO~12(volts)2/Hz 


(8.3) 


Although  this  noise  power  is  rather  formidable  by  itself,  it  is 
considerably  attenuated  at  16  KHz.  The  reason  is  that  this  resistor  and  the 
sensor  capacitor  act  as  an  RC  lowpass  filter  for  the  Nyquist  noise. 

The  effective  capacitance  to  ground  of  this  filter  is  that  of 
the  two  halves  of  the  sensor  in  parallel.  The  effective  lowpass 
cutoff  frequency  is  calculated  to  be  ^  *  1.4*102  Hz.7 


All  the  experimental  information  in  the  sensor  signal  is 

contained  in  a  small  bandwidth  about  the  reference  frequency  of 

16  KHz.  Synchronous  detection  shifts  this  bandwidth  nearly  tc  DC 

and  Nyquist  noise  outside  this  bandwidth  is  eliminated  by  the 

filtering  in  the  lock-in-amplifier.  Thus,  the  only  Nyquist  noise 

important  at  the  output  is  that  near  16  KHz.  At  this  frequency, 

the  lowpass  filter  . 

‘  has  an  attenuation  of 


1 

* 
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1  4x10  -3 

— — =■ —  =  1.4x10  in  amplitude.  The  noise  power 
10  _  .  •••*.•  - 
at  t*e  preamplifier  input  at  16KH,  is  reduced  by  the  square  of 

this  factor.  With  the  4  db  noise  figure  of  the  preamplifier  in 

the  present  circuit  «-aken  into  account,  the  effective  Nyquist 

noise  at  the  preamplifier  input  is 

P  (in)  =  8.7xl0'18  (volt)2/Hz  (8.4) 

v 


throughout  the  experimental  bandwidth.  At  the  output  of  the 
phase  sensitive  detector,  the  corresponding  noise  power  density 
in  the  region  of  flat  response  is 

p  .=  c2  P  (in)  *  8.7xlO'10(volt)2/Hz  (8.5) 

v  t  sens )  L  v 

Here  C. ,  the  DC  gain  of  the  preamplifier  and  phase  sensitive 
detector,  is  taken  as  lO4 ,  a  typical  value.  Comparison  of  the 
above  value  with  values  in  Table  3  shows  that  the  noise  level  at 
the  position  sensor  output  due  to  Nyquist  noise  is  at  least  an 
order  of  magnitude  below  that  due  to  Brownian  motion.  Thus,  Nyquist 
noise  presents  no  problems  in  practice. 

There  are,  however,  other  types  of  electronic  noise  in 

the  sensor  circuitry.  Noise  originating  in  the  mixer,  or  detector, 
of  the  phase  sensitive  detector  can  be  appreciable.  The  level  of 
this  noise  was  reduced  by  lock-in  amplifier  modifications  suggested 
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by  the  manufacturer.7  In  the  prefilters,  which  operate  at 
frequencies  near  DC,  1/f  noise  and  fluctuations  in  temperature 
and  supply  voltage  present  problems. 

Measurement  of  electronic  noise  in  all  the  instrument 
systems  is  accomplished  quite  simply.  The  drive  signal  to  the 
outer  plates  of  the  sensor  capacitor  is  removed  and  all  the 
electronic  circuits,  including  the  temperature  control  system,  are 
left  in  normal  operating  configuration.  Then  the  various  prefilter 
outputs  are  recorded  as  usual.  In  this  manner,  noise  from  all 
electronic  sources  is  observed,  while  Brownian  noise  and  ground 
noise  are  eliminated.  Root-mean-square  values  of  the  voltage 
fluctuations  at  the  prefilter  outputs  with  drive  removed  have 
been  computed  from  the  digitized  data.  They  are  as  follows: 


Normal  mode  prefilter 
(30  cph) 

V  (RMS)  2  5  mv 
out 


Seismic  prefiltcr  with  peak  gain  at  621 
Vout(RMS>  S  18  mv  (8.6) 


These  values,  of  course,  are  the  same  for  both  the  vertical  and 

horizontal  accelerometers,  as  the  electronic  systems  are  the  same. 

The  criterion  for  optimization  of  the  electronic  systems  is 

that  output  fluctuations  due  to  electronic  noise  from  all  sources 

be  smaller  than  those  due  to  Brownian  motion  in  the  accelerometer. 

If  this  is  achieved,  then  further  reduction  in  the  electronic  noise 

is  superfluous.  Comparison  of  results  8.6  with  Brownian  noise 

estimates  in  Table  4  shows  that  this  condition  is  nearly  achieved 

in  the  vertical  instrument.  Fluctuations  of  the  seismic  prefilter 

output  due  to  Brownian  motion  are  more  than  twice  as  large  as  those 

due  to  electronic  noise.  At  the  normal  mode  prefiltcr  output, 
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fluctuations  from  the  two  sources  are  nearly  equal.  In  the  horizontal 
instrument,  which  has  a  lower  Brownian  noise  level,  the  optimum  condi¬ 
tion  is  not  attained.  The  amplitude  of  electronics  noise  is  a  factor 
of  3  too  large  at  the  seismic  prefilter  output,  and  nearly  an  order 
of  magnitude  too  large  at  the  modal  prefilter  output.  Horizontal 
ground  noise,  however,  is  so  large  that  the  electronics  noise  level 
has  never  been  observed  with  drive  applied. 
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9*  Torsional  Mode  Data 

In  the  past,  the  torsional  free  oscillations  of  the  earth 
generated  by  large  earthquakes  have  been  measured  v/ith  both 
horizontal  accelerometers  and  strain  meters . 1 ‘Compared  to  the 
spheroidal  modes,  the  torsional  modes  of  the  earth  involve  larger 
fractional  shear  energies.  Because  shear  motion  is  more  dissipa¬ 
tive  than  compression,2  the  torsional  modes  have  lower  Q  .  Spher¬ 
oidal  fundamentals  typically  have  Q’s  near  300,  while  the  toroidal 
fundamentals  have  Q*s  of  about  150.  In  addition,  due  to  local 
tilting  of  the  ground,  horizontal  ground  noise  at  the  surface  in 
the  normal  mode  band  is  relatively  high. 

Because  of  the  higher  horizontal  noise  level  at  surface 
sites  and  the  generally  lower  Q  values,  torsional  oscillations 
are  observed  at  surface  sites  to  decay  into  the  noise  more  rapidly 
than  spheroidal  ones.  Thus,  observation  of  torsional  modes 
requires  large  initial  excitation  if  adequate  spectral  resolution 
for  separation  and  identification  of  modal  lines  is  to  be  obtained. 

Since  its  maiden  run  in  early  April,  1971,  the  horizontal 
accelerometer  has  been  operating  nearly  continuously  in  the  open 
loop  configuration.  For  most  of  this  time,  it  was  deployed  at  a 
surface  site  (vault  depth  10  ft.)  at  the  Elliott  Tield  Station 
(32. 88°N,  117. 1°W)  about  13  miles  northeast  of  San  Diego,  California. 
While  the  instrument  was  at  this  station,  no  particularly  large  cents 
occurred.  Such  is  the  capriciousness  of  nature,  however,  that  there 
were  two  major  events  during  the  few  days  required  to  transport  the 
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instrument  to  Wallace,  Idaho,  for  tests  in  the  Star  Mine  to 
determine  the  depth  dependence  of  earth  noise.  Thus,  no  observa¬ 
tions  of  strong  excitation  of  torsional  modes  were  obtained. 

To  illustrate  the  quality  of  data  available  from 
the  horizontal  accelerometer,  records  of  the  normal  mode  prefilter 
output  for  two  events  of  intermediate  magnitude  are  shown  in  Figures 
8  and  9.  Relevant  parameters  for  these  events  are  shown  in  Table  5. 


TABLE  5 


Parameter 

Origin  Time 
Origin  Date 
Epicenter 
Location 

Source-received  distance  A 
Incidence  Angle  a 
Magnitude 
Depth 


Adak  Earthquake 
(Figure  8) 

06:00:27.3  GMT 

May  2,  1971 

51 . 4°N ,  177. 2°W 

Adak  Island  Region 

46.8° 

93.1° 

Ms  =  7.1 
43  Km 


Chilean  Earthquake 
(Figure  9) 

21:00:40.9  GMT 

June  17,  1971 

25 . 5°S,  69 . 2°W 

Chile 

73.8° 

90.7° 
mfc  ~  6.3 
93  Km 


Except  for  the  incidence  angle  a  and  source-receiver  distance  A  , 
which  were  calculated  from  the  epicentral  locations,  the  information 
is  as  reported  in  the  National  Oceanographic  and  Atmospheric  Administi i 
tion's  Preliminary  Determinations  of  Epicenters.  The  incidence  angle 
o  is  the  smallest  angle  between  the  great  circle  path  from  the 
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FIGURE  8 


"°dC  Premt'r  data  <*«  loop  horizontal  accele,™- 
av  Elliott  Station  for  Adak  earthquake  of  May  2,  1971. 
Vertical  scale  applies  only  to  signals  in  1-30  cph  passband  of 
profilter,  start  time  of  record  is  06:12:30  GMT  May  7,  1971. 
Note  body  wave  arrivals  at  beginning  of  record  and  high  noise 

level  throughout.  Low  frequency  noise  is  probably  due  to 
local  ground  tilt 
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FIGURE  9 


Normal  mode  prefilter  data  from  open  loop  horizontal  accelerom¬ 
eter  at  Elliott  Station  for  Chilean  earthquake  of  June  17, 

1971.  Vertical  scale  applies  only  to  signals  in  1-30  cph  pass- 
band  of  prefilter.  Start  time  of  record  is  21:14:50  GMT  June  17, 
1971.  Note  large  body  wave  arrivals  at  beginning  of  record 
and  passage  of  Love  wave  packets  G1  and  G2  . 


epicenter  to  Elliott  St.tion  end  the  direction  of  maximum  sensitivity 

of  the  horizontal  instrument.  for  these  events,  the  horizontal  .cool 
cromctcr  was  oriented  HMS°E. 

The  f items  for  both  of  the  events  show  excitation  of  surface 
waves  which  last  for  several  hours.  In  neither  record,  however,  is 
significant  energy  observed  to  travel  completely  around  the  earl!,. 
Sine,  the  incidence  angle  o  is  close  to  90«  for  both  events,  the 
horizontal  notion  observed  is  due  to  transversely  polarized  Love  waves 
rather  than  longitudinally  polarized  Rayleigh  waves 

Thus  the  observed  disturbances  or.  interpreted  in  tcros  of  torsional 
free  oscillations  of  the  earth,  rather  than  spheroidal  oscillation:,, 
which  are  also  seen  on  horizontal  instruments.^ 

rieurc  10  shows  a  Fourier  power  spectiun  of  10  hours  of  the 
nodal  prefilter  record  for  the  Chilean  earthquake.  The  running  mean 
of  the  ambient  spectrum  is  plotted  for  comparison.  Below  lb  eph,  the 
signal-to-noisc  ratio  is  close  to  unity.  Above  this  fre,uency, 
definite  structure  appears  in  the  earthquake  spectrum  with  up  to  a 
20  db  signal-, o-noise  ratio.  Certain  lines  in  the  spectrum  are 
identified  as  fundamental  torsional  modes  of  free  oscillation  on  the 
basis  of  theoretical  eigenfreoueney  ealcul  tions.«  Although 
collective  motion  of  the  earth  as  a  whole  is  not  really  excited  by 
this  small  earthquake,  interpretation  of  the  disturbances  in  terms 
of  a  superposition  of  modes  is  still  possible.  The  relatively  high 

frequency  at  which  the  spectrum  emerg.es  from  the  noise  is  evidence 
of  the  rapid  decay  of  the  motion. 


312 


FIGURE  10 


Spectra  of  Chilean  earthquake  and  of  ambience  from  normal  mode 
prefilter  output  for  horizontal  accelerometer  at  Elliott  Station. 
Earthquake  spectrum  is  taken  from  10  hour  record  of  prefilter 
output,  of  which  section  shown  in  Eigu'  9  is  initial  part. 
Ambient  spectrum  is  derived  from  record  of  prefilter  output  of 
equal  length  beg'  .ning  06:24:59  GMT  April  30,  1971,  in  which 
there  was  no  vis  ble  earthquake.  Trace  shown  is  running  mean 

of  ambient  spectrum  over  1  eph  interval.  Boxcar  time  windows 
are  used  for  both  spectra.  Certain  toroidal  modes  of  earth 
are  tentatively  identified  from  theoretical  calculations  of 
eigenfrequcncies .  Eumpinesc  of  running  mean  of  ambient  spectrum 
is  due  to  high  horizontal  noise  level. 
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10.  Teleseismic  Detection  with  both  Horizontal  and  Vertical  Accel¬ 


erometers 

As  pointed  out  in  section  9,  horizontal  accelerometers  at 
surface  sites  are  subject  to  a  rather  high  noise  level  due  to  short- 
range  tilting  of  the  ground.  Since  the  amplitude  of  horizontal  dis¬ 
placement  in  Rayleigh  waves  is  only  about  2  to  3  times  that  of  verti¬ 
cal  displacement,  the  signal-to-noise  ratio  for  Rayleigh  waves  is 
higher  for  vertical  accelerometers.  Thus,  for  Rayleigh  wave 
detection  the  vertical  accelerometer  is  the  preferred  instrument. 
Nevertheless,  the  horizontal  instrument  is  quite  useful  for  seismic 
monitoring.  Due  to  their  horizontal  polarization,  Love  waves, 
which  provide  a  valuable  secondary  source  discriminant,  are  not 
observed  with  vertical  instruments.  In  addition,  due  to  the 
definite  horizontal  polarization  in  both  Rayleigh  and  Love  waves, 
the  horizontal  accelerometer  is  sensitive  to  the  direction  of  wave 
arrival.  Thus,  two  horizontal  instruments  can  be  used  to  estimate 
the  direction  of  a  source. 

In  early  April,  1971,  the  prototype  horizontal  quartz  fiber 
accelerometer  was  deployed  at  the  Elliott  Station  along  with  a 

vertical  accelerometer  The  two  instruments  were  placed  side-by- 

side  in  the  surface  vault  and  both  were  run  in  the  open-loop 
configuration.  For  comparison  wich  a  laser  strain  meter  at  the 
same  site,1 “♦  the  horizontal  was  oriented  N45°E.  Outputs  from  both 
instruments  were  recorded  digitally  on  magnetic  tape  at  a  sampling 
interval  of  1  second. 

Although  data  from  all  three  prefilterr.  were  continuously 
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recorded  for  both  instruments,  only  records  from  the  seismic 
prefilters  are  of  interest  here.  The  acceleration  response  of 
the  horizontal  instrument  at  the  seismic  prefilter  output  is 
plotted  relative  to  the  DC  open  loop  gain  in  Figure  7.  In  the 
seismic  band  below  the  resonance  frequencies  of  both  instruments, 
the  response  of  the  vertical  accelerometer  at  its  seismic  prefilter 
output  has  the  same  shape,  but  the  gain  level  is  different. 

The  DC  open  loop  gain  of  each  instrument  for  this  run  has 

been  carefully  measured.  For  the  vertical  accelerometer,  it  has 
been  doternined  by  tilt  calibration  and  by  a  fit  of  the  theoretical 
earth  tides  to  the  tidal  prefilter  output.  The  result  is 

GQ(vert)  =  (6.25  i  0.06)  x  106  volts/g  (10. 1] 

where  g  is  the  local  acceleration  of  gravity.  At  this  gain, 
peak-to-pcak  tidal  acceleration,  or  about  3  x  10_7g,  amounts  to  one 
tenth  of  the  full  scale  range  of  1  lOv  at  the  detector  output. 

For  the  horizontal,  the  DC  open-loop  gain  has  been  determined  by 
tilt  calibration  and  checked  against  design  calculations.  It  is 

^(hor)  -  (1.95  i  0.05)  x  10  volts/g  (10.2) 

Due  to  a  difference  in  the  attenuators  at  the  two  seismic  prefilter 
outputs  for  this  run,  0  db  for  the  horizontal  seismic  prefilter  in 
Figure  7  is  2.83  x  io  volts  per  g  of  acceleration  or  per  radian 
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of  tilt.  For  all  the  other  outputs  of  the  horizontal,  0  db  is 
1.95  x  10^  volts/g  as  above.  Thus,  at  its  seismic  prefilter  output, 
the  equivalent  acceleration  response  of  the  horizontal  is  a  factor 
of  2.2  less  than  that  for  the  vertical. 

Figure  11  shows  the  response  of  the  vertical  to  ground 

displacement  at  the  seismic  prefilter  output.  It  is  computed  using 
the  standard  seismic  approximation  in  which  only  kinematic  accelera¬ 
tion  is  considered.  Since  the  horizontal  instrument  measures  both 
kinematic  acceleration  and  tilt,  interpretation  in  terms  of  ground 
displacement  is  not  possible  for  it.  The  response  at  its  seismic 
prefilter  output  is  best  obtained  from  Figure  7. 

The  upper  two  traces  of  Figures  12  and  13  are 

computer  plots  made  from  simultaneous  recordings  of  the  seismic 
prefilter  output  for  both  instruments.  Trace  (a)  is  the  horizontal 
seismic  prefilter  output,  and  trace  (b)  the  vertical.  The  figures 
show  the  noise  levels  at  various  times  of  night,  as  well  as  seven 
distant  earthquakes.  These  two  examples  are  representative  of  a  much 

.larger  body  of  data  ,  part  of  which  has  been  reported  elsewhere.  5 

Careful  inspection  of  these  records  shows  that  the  noise 
from  the  horizontal  instrument  is  concentrated  at  frequencies 
slightly  below  those  of  seismic  interest.  In  fact,  most  of  the 
horizontal  noise  can  be  removed  from  the  records  by  digital  filter¬ 
ing.  The  lower  two  traces  in  Figures  12  and  13  show 

the  result  of  the  application  of  a  highpass  convolution  filter  to 
the  data  of  the  upper  two  traces.  This  highpass  filter,  whose 
response  is  shown  in  Figure  11,  has  a  -3  db  break  point  at  a 
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FIGURE  11 


Response  to  ground  displacement  at  seismic  prefiltcr  output  of 
open  loop  vertical  accelerometer  (scale  at  left),  and  response 
of  digital  highpass  filter  (scale  at  right).  Seismic  prefilter 
output  for  horizontal  instrument  is  not  readily  interpreted  in 
terms  of  ground  displacement.  Small  buinp  in  highpass  filter 
response  is  first  sideband.  Highpass  filter  response  is  essen¬ 
tially  unity  at  periods  shorter  than  45  seconds. 
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FIGURE  12 

Computer  plots  of  horizontal  and  vertical  seismic  prefilter  outputs 
for  Event  4  of  Table  6.  Traces  (a)  and  (b)  are  raw  outputs 
of  horizontal  and  vertical  seismic  prefilters,  respectively. 

Traces  (c)  and  (d)  are  the  results  of  application  of  the  digital 
filter  of  Figure  11  to  traces  (a)  and  (b). 
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FIGURE  13 

Computer  plots  of  horizontal  and  vertical  seismic  prefilter  outputs 
for  Event  5  of  Table  6.  Traces  (a)  through  (d)  are  as  in  Figure  12. 
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period  of  45  seconds  and  is  down  20  db  at  50  seconds.  After 
digital  filtering,  the  peak-to-peak  horizontal  noise  (trace  (c)) 
is  reduced  by  nearly  an  order  of  magnitude,  yet  the  seismic  signal 
unaffected  by  the  filtering.  The  signal-to-noise  ratio  for  the 
vertical  record  is  also  improved  by  digital  filtering,  as  shown  in 


trace  (d). 

The  surface  site  at  Elliott  Station  is  subject  to  consi¬ 
derable  acoustic  and  high  frequency  ground  noise  from  a  nearby  air 
base  and  a  highway.  However,  noise  in  the  seismic  band  is  rela¬ 
tively  low,  at  least  for  the  vertical.  The  horizontal  noise 

is  much  higher  and  varies  with  weather  and  the  time  of  day. 
If  the  output  of  the  horizontal  is  interpreted  in  terms  of  accelera¬ 


tion,  and  if  the  difference  in  DC  gain  between  the  two  instruments 
is  taken  into  account,  the  vertical  and  horizontal  noise  levels 
can  be  compared.  At  night,  the  horizontal  noise  is  a  factor  of 
6  to  15  greater  in  equivalent  acceleration  than  the  vertical. 

During  midday,  this  factor  is  augmented  by  as  much  as  another  factor 
of  three,  especially  during  windy  and  variably  cloudy  days.  When 


the  weather  is  calm  and  uniformly  overcast,  as 
daytime  horizontal  noise  is  as  low  as  it  night. 


on  foggy  days,  the 
This  variation  in 


the  horizontal  noise  level  is  due  to  the  fact  that  the  noise 
derives  from  short-range  pressure  and  temperature  fluctuations .9 , 18 


Noise  originating  in  the  instrument  electronics  has  been 
measured  by  removing  the  reference  signal  as  explained  in  section 
8.  The  seismic  prefilter  output  for  either  instrument  with  no 
reference  makes  a  straight  line  on  the  scales  of  either  of  Figures  12 
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and  13. 


The  effects  of  Brownian  motion  have  been 
calculated  for  each  instrument  for  this  run  in  the  manner  of  section 
9.  A  the  outputs  of  the  seismic  prefilters,  the  voltage  fluctua¬ 
tions  due  to  Brownian  noise  have  an  amplitude  of  38  mv  for  the 
vertical  and  3.5  mv  for  the  horizontal.*  As  Figures  12  and  13 
show,  this  is  at  least  an  order  of  magnitude  down  from  the  noise 
observed.  Coherence  measurements  with  a  second  horizontal 
will  be  made  to  determine  total  instrumental  noise  levels.  At 
present,  the  preliminary  noise  measurements  and  calculations  above, 
as  well  as  the  variation  of  the  horizontal  noise  with  environmental 
conditions,  lead  to  the  conclusion  that  the  observed  noise  is 
predominantly  real  earth  motion. 

Table  6  gives  the  identification  and  parameters  of 
the  events,  which  occurred  on  April  14  and 

April  17 ,  1971.  The  angle  a  in  the  table  is  the  smallest  angle 
between  the  great  circle  path  connecting  the  epicenter  to  the 
station  and  the  line  of  maximum  sensitivity  of  the  horizontal. 

When  this  angle  is  near  90°,  only  transversely  polarized  Love 
waves  are  observed  by  the  horizontal.  When  it  is  close  to  zero, 
only  Rayleigh  waves,  which  are  longitudinally  polarized,  are  seen. 

*For  the  horizontal,  this  number  is  slightly  different  from  that 
given  in  Table  4  because  different  prefilter  attenuation  was  used 
for  these  observations. 
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The  origin  times,  locations,  epicenters,  body  wave  magnitudes  and 
depths  of  the  events  are  as  given  in  the  L.A.S.A.  or  N.O.A.A.  P.D.E. 
reports.  The  source-receiver  distance  A  and  the  angle  a are 
computed  from  the  reported  epicenters. 

Both  the  surface  wave  magnitudes  M  and  M  (30°)  are  calculated  from 

s  s 

the  vertical  records  in  Figures  12  and  13,  using  the  calibration 
implicit  in  Figure  11.  The  method  is  described  in  detail  elsewhere  ,7» 1 14 
The  surface  wave  magnitude  Ms(30°)  is  the  magnitude  computed  by 
setting  the  source- receiver  distance  arbitrarily  to  30°.  It  is 
useful  in  estimating  the  effectiveness  of  instrumentation  for  seismic 
monitoring  of  underground  nuclear  explosions.  Detection  capability 
of  the  present  instrumentation  is  examined  in  detail  elsewhere.5*  7»  1** 
Figure  12  shows  Event  4,  whose  seismic  waves  arrived  at  an 
angle  of  78°  to  the  horizontal's  axis  of  maximum  sensitifity.  Strong 
Love  waves  are  observed,  so  that  the  signal--to-noise  ratio  is 

appreciably  larger  in  the  horizontal  records.  Figure  13  shows 
Event  5,  whose  seismic  waves  arrived  at  an  angle  of  39°  to 
the  horizontal's  axis  of  maximum  sensitivity.  No  Love  waves  were 
observed  for  this  event,  and  the  signal-to-noise  ratio  was  6  times 
larger  on  the  vertical.  There  is  an  apparent  decrease  of  the 
horizontal-to-vertical  polarization  ratio  of  the  Rayleigh  waves 
with  increasing  period  in  this  record. 

Although  the  new  wide- band  horizontal  accelerometer  is 
noisy  at  this  surface  site,  the  noise  is  not  in  the  frequency 
region  of  seismic  interest.  It  can  be  removed  by  digital  filtering. 

The  digitally  filtered  horizontal  records  then  become  a  powerful 
adjunct  to  the  vertical  records  for  purposes  of  seismic  detection. 

Both  the  theoretical  power  of  Love-waves  as  a  source  discriminant 
and  an  augmented  signal-to-noisc  ratio  for  events  producing  Love 
waves  are  now  available.  84 
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